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Motivated by a variety of theories that predict new effects, we tested the gravitational 1/r 2 law at 
separations between 10.77 mm and 137 using two different 10-fold azimuthally symmetric torsion 
pendulums and rotating 10-fold symmetric attractors. Our work improves upon other experiments 
by up to a factor of about 100. We found no deviation from Newtonian physics at the 95% confidence 
level and interpret these results as constraints on extensions of the Standard Model that predict 
Yukawa or power-law forces. We set a constraint on the largest single extra dimension (assuming 
toroidal compactification and that one extra dimension is significantly larger than all the others) of 
R, < 160 /jm, and on two equal-sized large extra dimensions of R* < 130 /xm. Yukawa interactions 
with \a\ > 1 are ruled out at 95% confidence for A > 197 /jm. Extra-dimensions scenarios stabilized 
by radions are restricted to unification masses M* > 3.0 TeV/c 2 , regardless of the number of large 
extra dimensions. We also provide new constraints on power-law potentials V(r) oc r~ k with k 
between 2 and 5 and on the 75 couplings of pseudoscalars with m < 10 meV/c 2 . 

PACS numbers: 04.80.Cc,04.80.-y 



I. INTRODUCTION 
A. Background 

Until a few years ago, it was widely assumed that the 
Newtonian Inverse Square Law (ISL) should be valid for 
length scales from infinity to roughly the Planck length 
Rp = y/ Gh/c 3 = 1.6 x 10" 35 m, at which scale quan- 
tum effects must become important. After all, the usual 
argument went, the exponent 2 in the force law simply 
reflects the fact that we live in a 3-dimensional world. 
A wide variety of recent theoretical speculations, moti- 
vated in part by string-theory considerations, have raised 
the possibility that fundamentally new phenomena could 
occur at length scales below 1 mm. Many of these specu- 
lations are driven by the two so-called hierarchy problems 
of gravity: 

• The gauge hierarchy problem. Gravity is extraor- 
dinarily weak compared to the other fundamental 
forces. The Planck mass Mp = yJhc/G = 1.2 x 
10 16 TeV/c 2 is huge compared to the electroweak 
scale Mew ~ 1 TeV/c 2 . It has been argued[l] 
that the true Planck mass, M#, could be as low as 
1 TeV/c 2 if some of the "extra" space dimensions 
demanded by string-theory are "large" compared 
to the Planck length. It is possible that the size 
of some of the "large extra dimensions" could be 
large enough to alter the gravitational Gauss law, 
so that gravity would become anomalously strong 
in an experimentally accessible regime [1]. 



small compared to the predictions of quantum me- 
chanics. The gravitating energy density p vac ~ 
0.7p c , inferred from a wide variety of astrophysi- 
cal observations [2], is at least 10 60 times smaller 
than the predicted zero-point energy for a cutoff 
of Mp. The observed energy density corresponds 
to a length scale i? V ac = \/hc/ Pv&c ~ 0.1 mm and 
an energy of \J (hc) 3 p vac w 2 meV that may have 
fundamental significance [3]. It has been suggested 
that the apparent inability of gravity to "see" the 
vacuum energy could be explained if the effective 
theory of gravity had a cutoff of ~ 1 meV[4, 5], so 
that gravity would effectively "shut off" at length 
scales less than i? V ac- 

Experimental tests of the gravitational ISL also probe 
certain speculations about non- gravitational physics. 
The standard model of particle physics cannot be com- 
plete and many ideas for extending it predict very-low- 
mass scalar or vector bosons that could produce short- 
range exchange forces that would appear as violations 
of the ISL (see, for example, an extensive summary in 
Rcf. [6]). 

The desire to test a basic law in a previously inaccessi- 
ble, but very interesting, regime motivated the work we 
report here. Some of the work we report in this paper has 
already appeared in Letter form[7]. This paper includes 
additional experimental work and an improved analysis; 
it supercedes Ref. [7]. 



B. Parameterizations 



• The cosmological constant problem. The observed 
gravitating vacuum-energy density is vanishingly 
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It is now customary to interpret experimental tests of 
the ISL as setting bounds on a possible Yukawa addition 
to the familiar 1/r Newtonian potential 
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where a is a dimensionless strength parameter and A a 
length scale. The Yukawa potential is, of course, the 
static limit of the interaction from exchange of a bo- 
son of mass m = h/(cX) in which case a is proportional 
to the squared product of the appropriate coupling con- 
stants. This Yukawa form is obviously appropriate for 
the boson-exchange forces mentioned above. It is also a 
good approximation to the effects one expects from large 
extra dimensions until the separation of the interacting 
bodies becomes comparable to or smaller than the size of 
the large extra dimensions [8]. 

The ISL can also be violated by power-law potentials, 
which we parameterize as 
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these can occur in processes where two massless parti- 
cles are exchanged. We therefore interpret our results as 
exclusion plots in the a — A parameter space and as con- 
straints on power-law parameters (3k for k = 2, 3, 4 and 
5. 



C. Previous Experimental Limits 

The Yukawa constraints on ISL violation for length 
scales A < 1 cm, at the time this work was begun in 
late 1999, are shown in Fig. 1. All the constraints in 
this Figure were obtained from torsion-balance experi- 
ments: a null test using cylindrical geometry from the 
Irvine group [9] set the most restrictive limits for length 
scales larger than about 2.5 mm, an absolute comparison 
in Moscow of the force between spheres at two differ- 
ence distances set the best limit for 80 /im < A < 2.5 
mm[10] and a measurement of the Casimir force between 
a sphere and a plane[ll, 12] set the tightest constraint 
for 4 /im < A < 80 /jm. 



II. THEORETICAL SPECULATIONS 

A surprisingly rich array of theoretical speculations 
suggests that fundamentally new physics could show up 
at length scales between 1 mm and 100 /im. This new 
physics is clearly probed directly by ISL tests, but may 
also be revealed in collider experiments at the very high- 
est energies or in astrophysical phenomena such as neu- 
trino production in supernova explosions. The relative 
sensitivity of these very different probes depends strongly 
on the basic assumptions of specific models. We give here 
only a brief summary of these considerations and refer the 
reader to two recent reviews[6, 13] for more comprehen- 
sive discussions. 
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FIG. 1: [color online] 95% confidence upper limits on a 
short-range Yukawa violation of the gravitational ISL as of 
1999. The region excluded by previous torsion-balance exper- 
iments [9-12] lies above the heavy lines labeled Lamoreaux, 
Moscow, and Irvine, respectively. The constraints are taken 
from Ref. [14]. The dilaton and moduli predictions are taken 
from Refs. [15, 16], respectively. The extra dimensions and 
radion predictions are based on Eqs. 3, 4 and 6, assuming 
n = 2 and A'/, = 1 TeV. 



A. Large Extra Dimensions 



The discovery that string or M-theory necessarily con- 
tains objects called branes, where the number of spatial 
dimensions of the brane can be anywhere from to 9, has 
opened new possibilities for explaining the gauge hierar- 
chy problem. Arkani-Hamed, Dimopoulos and Dvali[l] 
proposed that the "true" Planck mass could be as low as 
1 TeV; in their scenario all the Standard Model particles 
are represented by open strings whose ends are fixed to a 
3+1 dimensional brane, while the graviton which is rep- 
resented by a closed loop of string is free to propagate 
in the 10+1 dimensional "bulk". In their picture gravity 
appears weak not because it is inherently weak, but be- 
cause most of its strength has leaked off into the extra 
(curled up) dimensions of the bulk which are inaccessi- 
ble to physicists living on the brane. This notion can 
be tested very directly by probing the gravitational ISL 
at length scales comparable to the radius of the curled 
up dimensions i?* which, it is suggested, could be much 
larger than the Planck length. If the effective number of 
large extra dimensions is n, the gravitational Gauss-law 
force will smoothly change from a 1/r 2 form for r >> i?,* 
to a \/r 2+n form for r « i?*. For the experimentally 
relevant case where r > the gravitational potential 
is expected [6, 17], assuming compactification on a torus, 
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to have the form of Eq. 1 with 

A = i?* and a = 8n/3 , 



(3) 



so that gravity is predicted to get much stronger when 
probed at distances comparable to i?*. If all n "large" 
extra dimensions have the same size i?*, conservation of 
flux gives a relation between the "true" Planck mass M* , 
i?*, n and the Newtonian constant G. For compactifica- 
tion on a torus 
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Although the case with n = 1, ii* = 3 x 10 12 m is ruled 
out by solar-system observations, at the time this work 
was begun the case with n = 2, R* = 0.3 mm, as well as 
scenarios with n > 2, were consistent with ISL data. 



B. New Boson- Exchange Forces 

The character of the force produced by boson exchange 
depends on the spin and parity of the exchanged parti- 
cle. Natural-parity (0 + or 1~) bosons give forces be- 
tween unpolarized bodies, while unnatural-parity (0 _ or 
1 + ) bosons give purely spin-dependent forces. The forces 
from boson exchange violate the Weak Equivalence Prin- 
ciple and can be distinguished from forces due to extra 
dimensions by Equivalence Principle tests. 



1. Radion-mediated Forces 

In string theories, the geometry of spacetime is ex- 
pected to be dynamical with the radii of new dimen- 
sions fluctuating independently at each point in our 4- 
dimcnsional spacetime. In an effective low-energy theory, 
the volume of the extra dimensions must be stabilized 
by radions, low-mass spin-0 fields with gravitational- 
strength couplings that determine the radius of the new 
dimensions. Radion exchange will produce a force with 
a strength [6, 17] 



n + 2 



and range of order 



A 



cGMi 



2.4 



1 TcV 



M»c 2 



(6) 



(7) 



In many cases the radion-mediated force is the longest- 
range effect of new dimensions[18], and one which, unlike 
the direct effect (Eq. 5), does not diminish as the number 
of new dimensions increases. 



2. Exchange of String- Theory Scalars 

All known acceptable ground states of string theory are 
supersymmetric and contain large numbers of extremely 
weakly coupled, massless scalar fields called moduli. The 
expectation values of the moduli fields set the parame- 
ters of the effective theory. The moduli necessarily cou- 
ple weakly to the supersymmetry-breaking sector, and 
for a low supersymmetry-breaking scale are expected to 
acquire very small masses[16] . Moduli couplings are com- 
putable in any given vacuum, so that ISL tests are the 
best way to search for these particles. The best under- 
stood of these predicted scalar particles is the dilaton 
which determines the strength of gauge couplings. Its 
couplings to ordinary matter arc nearly free of QCD 
uncertainties[15] so that its discovery could provide a 
smoking gun for string theory. 



3. Axion- Exchange Force 

The force due to exchange of axions, pseudoscalar par- 
ticles invented to explain the very small upper limit on 
©QCD , is interesting. The first-order force from exchange 
of unnatural-parity bosons is purely spin-dependent and 
vanishes between unpolarized objects. However, if the 
axion does not drive Oqcd all the way to zero, but merely 
to some very small value, the axion acquires a small, 
CP-violating, scalar admixture. This, in turn, leads to 
a spin-independent Yukawa potential between nucleons 
with[19] 
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where m a is the axion mass which is constrained to lie 
between ss 1 /zeV/ c 2 and « 10 meV/c 2 [20] , and f n arc 
the pion mass and decay constant respectively, m u and 
rrid are the current masses of the up and down quarks re- 
spectively, and u is the atomic mass unit. The numerical 
value in Eq. 8 assumes m u /md = 0.5. The experimental 
upper limits on the neutron and mercury-atom electric 
dipole moments imply that Oqcd < 6 x 10~ 10 [21] and 
©QCD < 1-5 x 10~ 10 [22], respectively. 



4- Multi-Particle Exchange Forces 

Power-law potentials such as Eq. 2 can arise from the 
2nd-order process in which 2 massless particles are ex- 
changed simultaneously. Such processes are particularly 
interesting when the exchanged particles are unnatural- 
parity bosons (for which the lst-order force vanishes 
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when averaged over unpolarized test bodies) or fermions 
(for which the lst-order process is forbidden). The great 
sensitivity of "gravitational" experiments using unpolar- 
ized test bodies leads to useful constraints on the cou- 
plings in these cases. 

Potentials with k = 2 or k = 3 may be generated by the 
simultaneous exchange of two massless scalar [23] or 75- 
coupled massless pseudoscalar[24] particles, respectively. 
The exchange of two massive pseudoscalar bosons having 
75 couplings to fermions a and b with coupling constants 
g a and gb, produces a long-range, spin-independent ex- 
change potential[25, 26], 
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where K\ is a modified Bessel function, A = h/(mc), and 
M and m refer to the fermion and pseudoscalar masses 
respectively. In the limit where the pseudoscalar boson 
is massless Eq. 10 becomes 
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FIG. 2: [color online] Torsion pendulums and rotating attrac- 
tors used in Experiments I (left) and II (right). The active 
components are shaded. For clarity, we show an unrealis- 
tically large vertical separation between the pendulums and 
attractors, and omit the conducting membranes and attractor 
drive mechanisms. 



Potentials with k — 5 may be produced by the simultane- 
ous exchange of two massless pseudoscalar particles with 
couplings such as axions or Goldstone bosons[24], 
or by a massless neutrino-antineutrino pair [27]. 



III. APPARATUS 



A. General Principles 



C. Vacuum-Energy Scenarios 

The cosmological constant problem and theoretical at- 
tempts to solve it have been extensively reviewed (see for 
example Refs. [2, 28]). These efforts to solve the problem 
fall into 2 broad classes: attempts to find a mechanism for 
radically reducing the quantum mechanical prediction for 
the vacuum energy density, or attempts to find a mecha- 
nism for reducing the gravitational coupling to the stan- 
dard vacuum energy. The latter attempts are particularly 
interesting from the standpoint of this work. Beane[3] 
argued that in any local effective quantum field theory, 
naturalness implies new gravitational physics at length 
scales of about a millimeter that would cut off shorter dis- 
tance contributions to the vacuum energy. Sundrum[4] 
proposed that the graviton is a "fat" object with a size 
of about i? V ac and has been exploring how this might re- 
duce its coupling to the vacuum energy [5], although it 
is not yet clear how self-consistent this is. This scenario 
makes a definite prediction that gravity "shuts off" at 
length scales below about 100 fim. In the framework of a 
Yukawa ISL- violation, this corresponds to a = — 1 with 
A ^ 0.1 mm. 



To test the ISL at smaller length scales than had been 
studied before, we developed a new torsion-balance in- 
strument, shown schematically in Fig. 2, that used pla- 
nar test bodies rather than cylindrical[9] or spherical[10] 
bodies that had been employed previously. Our test bod- 
ies were the "missing masses" of holes bored into cylin- 
drically symmetrical plates. In each of the two first- 
generation experiments reported here, the active compo- 
nent of the torsion pendulum was a thin ring containing 
10 cylindrical holes equally spaced around the azimuth, 
and the pendulum was suspended above a uniformly ro- 
tating, circular attractor disk containing 10 similar holes. 
In the absence of the holes, the disk's gravity would pull 
directly down on the ring and could not twist it. But 
because of the holes, the ring was twisted by a torque 
N((f>) — —dV((f))/d(f) where V(4>) is the potential energy 
of the ring in the field of the disk when the disk's holes 
are displaced by an angle <f> with respect to those in the 
pendulum. This torque oscillated 10 times for every rev- 
olution of the disk. V(4>) was not a simple sinusoidal 
function of <f> so that rotating the attractor at frequency 
u> produced torques at frequencies of lOw and its integer 
multiples. The IOlj, 2Quj, and 30w torques, N w , N 2 o, 
N30, were measured as functions of the vertical separa- 
tion between the bottom of the pendulum and the top of 
the attractor (higher harmonic twists were greatly atten- 
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uated by the inertia of the pendulum and did not provide 
useful signals). By placing the signals at high multiples 
of the rotation frequency, u, we reduced many potential 
systematic errors. We minimized electrostatic interac- 
tions between the attractor and pendulum by interpos- 
ing a stiff conducting membrane between the attractor 
and the pendulum and surrounding the pendulum with 
an almost complete Faraday cage. 

The experiments were turned into approximate null 
measurements by attaching a second, thicker, disk to the 
bottom of each attractor. This disk also had 10 equally- 
spaced holes bored into it, but the holes were rotated by 
18 degrees compared to those in the upper disk. The di- 
mensions of these thicker and larger-diameter holes were 
chosen so that the 10w Newtonian torque on the pendu- 
lum from the upper attractor holes was essentially can- 
celed by the lOw Newtonian torque from the lower holes. 
On the other hand, torques from a short-range interac- 
tion with a length-scale less than the thickness of the 
upper attractor disk could not be canceled because the 
lower attractor was too far from the pendulum ring. The 
20uj and 30w Newtonian torques were reduced much less 
substantially than the Newtonian lOw torque. 

Data were taken at separations ranging from s = 10.77 
mm to s = 137 fim, where s is the distance from the top 
of the attractor to the bottom of the pendulum. The 
signatures distinguishing conventional gravity from new 
short-range physics were 

1. a characteristic shape of the Newtonian lOw torque 
N^. The cancelation of was a strong function 
of s. In Experiment I, N^ was exactly canceled 
at s i=a 2 mm, undercanceled for s < 2 mm and 
overcancelcd for s > 2 mm. On the other hand, 
for ranges of interest the Yukawa torque was a 
monotonically decreasing function of s. As a result, 
the exact location of the zero-crossing was very sen- 
sitive to any violation of the ISL. In Experiment II, 
the cancelation occured at s w 3.3 mm. 

2. a relatively high harmonic content of the Newto- 
nian torque. N§) and N^q were comparable to Yfj 
because N^ was highly canceled while N§) and N^ 
torques were not. On the other hand, N^q and 
will be much less important than N^ - 

The predicted Newtonian, Yukawa and power-law 
torques are shown as functions of s in Fig. 3. 

We inferred the harmonic components of the torque 
from the pendulum twist which we measured by reflect- 
ing a laser beam from a mirror attached to the pendu- 
lum body. These measured torques were then compared 
to calculations of the expected Newtonian and possible 
Yukawa and power-law effects, allowing for uncertain- 
ties in key experimental parameters, such as the actual 
masses of the holes. Agreement of the best-fit values of 
these parameters with values determined independently 
added to our confidence in the results. 
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FIG. 3: Predicted torques for Experiment I as a function 
of separation s. Solid curves are Newtonian gravity; dotted 
curves are for a Yukawa interaction with a — 3, A = 0.25 mm; 
dashed curves are for a power-law interaction with /?3 = 0.05. 



TABLE I: Properties of the pendulum holes. All quantities 
are averages over the 10 holes. The mass refers to the missing 
mass of a single hole. 



Property 


Experiment I 


Experiment II 


hp (mm) 


2.002 ± 0.002 


2.979 ± 0.002 


a p (mm) 


4.7725 ± 0.0015 


3.1873 ±0.0015 


Mp/10 (g) 


0.3972 ± 0.0004 


0.266177 ± 0.000014 


r p (mm) 


27.665 ± 0.005 


26.667 ± 0.005 



B. Torsion Pendulum 

The torsion pendulums used in the two experiments, 
shown in Fig. 2, were quite similar. The active compo- 
nents of the pendulums were rings made from 7075 alu- 
minum with 10 cylindrical holes evenly spaced about the 
azimuth. The relevant properties of the rings — the hole 
radius a p and height h v , the distance r p from the hole 
centers to the pendulum's symmetry axis, and the "miss- 
ing mass" M p (the mass removed by drilling the holes) 
are shown in Table I. The mass measurements are dis- 
cussed in Section VI C 1. All quoted masses include small 
corrections for buoyancy in air and for the thin layer 
of gold that covered the entire surface. This gold layer 
provided an uniform electrostatic surface that minimized 



TABLE II: Properties of the Experiment I attractor holes. All 
dimensions are averages over the ten holes; r a is the distance 
of the hole centers from the rotation axis, and the masses are 
given for a single hole. 

Property Upper Disk Lower Disk 

K (mm) 1.847 ± 0.005 7.828 ± 0.005 

a a (mm) 4.7690 ± 0.0015 6.3449 ± 0.0015 

M a /10 (g) 1.1770 ± 0.0004 8.8666 ± 0.0019 

r a (mm) 27.655 ± 0.005 27.664 ± 0.005 
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variations in the surface electrical potential; gold does 
not oxidize and has unusually small differences in the 
work functions of its different crystal faces (see Ref. [6]). 
The entire pendulums, including the mirrors, as well as 
the tops of the upper attractor plates were gold coated 
for this same reason. In addition, the pendulums were 
nested inside a gold-coated copper shroud that contained 
small holes for the torsion fiber and the laser beam. 

The pendulum rings were mounted on gold-coated alu- 
minum frames that held the mirrors used by the optical 
system which measured the pendulum twist. The frame 
for Experiment I had two-fold azimuthal symmetry with 
two mirrors. In Experiment II we reduced the gravity- 
gradient noise by adopting a four-fold symmetric frame 
with four mirrors. The upper portions of the frames con- 
tained screws for leveling the pendulum rings to the hori- 
zontal. In Experiment I, these screws translated the fiber 
attachment point relative to the rest of the pendulum 
body, while in Experiment II, the screws adjusted the 
pendulum's center-of-mass so that it lay directly under 
the fiber. The total mass suspended from the fiber was 
28.2527±0.0002 g in Experiment I and 39.6463±0.0002 g 
in Experiment II. 



C. Rotating Attractor 

The attractors for Experiments I and II, shown in 
Fig. 4, were pairs of flat, concentric, high-purity cop- 
per disks. The upper disks were thinner and contained 
holes which, when properly rotated, aligned with those on 
the pendulum. The torque produced by these "in-phase" 
holes on the pendulum was opposed by the 10 "out-of- 
phase" holes in the lower disks that were rotated by 7r/10 
relative to those in the upper plate (the upper attractor 
disk in Experiment II also contained out-of-phase holes). 
The disks for both experiments were lapped to a flatness 
of ±0.003 mm and the two sides of each disk were parallel 
to better than 30 /zrad. The top surface of the attractor 
assembly was coated with gold. Tables II and III show 
the relevant properties of the disks. 



FIG. 4: [color online] Top and side views of the attractor 
disks. The disks from Experiment I are on the left. The solid 
circles represent the in-phase upper-plate holes. The open 
circles show the out-of-phase holes in the upper and lower 
plates. The three small holes held pins that aligned the upper 
and lower disks. 



TABLE III: Properties of the Experiment II attractor holes. 
Uncertainties in the least significant figures are shown in 
parentheses. 



Property 


Upper 


Disk 


Lower Disk 




In-phase 


Out-of Phase 


Out-of-Phase 


h a (mm) 


3.005(2) 


3.005(2) 


6.998(2) 


a a (mm) 


3.1930(15) 


4.7720(15) 


7.9693(15) 


M a /10 (g) 


0.86228(4) 


1.92592(5) 


12.50984(6) 


r a (mm) 


26.685(5) 


38.697(5) 


36.666(5) 



D. Electrostatic Shield 

The main component of the electrostatic shield was 
a tightly stretched 20-/im-thick, 11.4-cm-diameter BcCu 
foil clamped between two lapped aluminum rings. The 
foil and rings were assembled in a liquid nitrogen bath. 
As the unit warmed to room temperature, the differential 
thermal expansion between the BeCu foil and the alu- 
minum rings created the high tension necessary to hold 
the foil flat. Inspection with a microscope showed that 
the membrane was flat to within ±15 /im. The tension 
was sufficient to give the foil a fundamental (m=0, n~l) 
drum-head frequency of /o,i ~850 Hz. We found this res- 
onance by raising the frequency of an acoustical source 
until small particles placed on the membrane began to 
jiggle- 

E. Pendulum Suspension System 

The pendulum suspension system consisted of a main 
torsion fiber hung from an eddy-current damper that, 
in turn, hung from a short "pre-hanger" fiber that was 
attached to a computer-controlled x-y-z-9 stage. The 
system was electrically insulated so that we could mea- 
sure the pendulum-electrostatic shield capacitance (see 
Sec. IV B 2). The pendulum was grounded externally to 
the electrostatic shield in normal operation. 
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Torsion Fiber 



Attractor Rotation Drive 



The main torsion fiber was an 80-cm-long, 20-/xm- 
diameter, gold-coated tungsten wire. The fiber was at- 
tached to the pendulum frame with a crimp tube embed- 
ded in a copper screw. The fiber could support a 100 g 
load, so the pendulums, which had masses less than 40 g, 
did not stress the fiber near its elastic limit. The fiber's 
torsional spring constant, k, may be determined from 
the relationship n = £1 2 /, where ft is the pendulum's 
resonant angular frequency and / is its rotational inertia 
which we computed from a detailed model of the pendu- 
lum based on the measured dimensions and masses of its 
components. 

In Experiment I, fl was 15.720±0.002 mrad/s and we 
calculated I = 128 ± 1 g-cm 2 , yielding k = 0.0317 ± 
0.0002 erg/rad. In Experiment II we found ft = 13.454 ± 
0.002 mrad/s and computed / = 175 ± 1 g-cm 2 , yield- 
ing k — 0.0317 ± 0.0002 erg/rad. As expected, the in- 
ferred k's were essentially identical because the same tor- 
sion fiber was used in both experiments. The fiber hung 
from a passive damper consisting of an aluminum shaft 
and disk suspended between two ceramic ring magnets. 
The damper's axial symmetry minimized any damp- 
ing of the torsional motion (which by the fluctuation- 
dissipation theorem would have contributed to our noise) , 
but the damper "killed" the swing, wobble and guitar- 
string modes of the pendulum; the typical exponential 
damping time of such oscillations being approximately 
1000 s. The top of the damper shaft, in turn, hung from 
a 3-cm-long, 150-/im-diameter "pre-hanger" fiber. This 
torsionally stiffer prehanger acted as a gimbal, reducing 
any tilt-twist coupling to the torsion fiber. 

The fiber was prepared by heating it under load in vac- 
uum to ~ 70°C for 12-24 hours. After a day of cooling, 
the fiber's drift rate-the slow change in the fiber's equi- 
librium twist angle-was typically reduced to < 1 /irad/h. 



2. Fiber Positioning Mechanism 

The pendulum was positioned relative to the attractor 
by hanging the suspension fiber from an x-y-z-6 stage 
located outside the vacuum chamber at the top of the 
instrument. The shaft holding the pendulum suspension 
system was fed through the vacuum vessel via a rotary O- 
ring seal, and a welded bellows allowed for horizontal and 
vertical translation. The full range of vertical translation 
was 25 mm, while the horizontal travel was limited to 8 
mm. 

Three electronically controlled actuators (Newport 
Corporation model CMA-25PP with an ESP-300 con- 
troller) drove the x, y, and z stages. The rotation stage 
(Newport Corporation model URM80APE) controlled 
the equilibrium angle of the torsion fiber and was also 
used to remove torsional energy from the pendulum be- 
fore starting a run. 



The attractor was held in a copper cup with three- 
fold azimuthal symmetry. The symmetry was designed 
to allow for gravitational centering (see, for example, 
Ref. [29]), but was ultimately not used as a more ef- 
fective method was found (see Sect. IV Bl). The cup 
was mounted on an electrically insulated and largely 
non-magnetic bearing consisting of an 8 mm-diameter 
stainless-steel shaft held by two bearings with ceramic 
balls riding in Stcllite races. The entire bearing assem- 
bly was mounted on 3 leveling screws that were used to 
make the attractor precisely parallel to the electrostatic 
shield as discussed in Section IV A 2. 

The attractor drive motor (made by SAIA-Burgess) 
was a 48-pole stepper with a 5000:1 reduction gear giving 
240,000 steps per attractor revolution. Drive pulses were 
provided by a Stanford Research Systems DS345 Func- 
tion generator and the attractor angle was determined 
by counting the drive pulses. The absolute angle was ob- 
tained from a LED /photodiode pair that was interrupted 
once per revolution by a tab on the attractor drive mech- 
anism, producing an index pulse that was recorded by our 
data-taking computer. The index was constant to within 
0.1 degrees throughout a run (~10 attractor revolutions). 
The motor was mounted outside the vacuum vessel and 
drove the attractor through an O-ring seal. 



G. Autocollimator 

The pendulum twist was measured by the autocollima- 
tor shown in Figure 5. Collimated light from a 20 mW, 
780 nm, actively stabilized laser diode (Sony SLD301), 
placed at the focus of a 30 cm lens, was reflected by a 
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FIG. 5: [color online] Schematic diagram of the optical system 
that measured the pendulum twist. 



8 



polarizing beam splitter before passing through the lens 
to produce a 6-mm-diamctcr parallel beam. The beam 
entered the vacuum vessel through an anti-reflective win- 
dow, was reflected twice by two fixed, flat mirrors, and 
then struck a flat mirror on the pendulum. The light re- 
flected off the pendulum struck a fixed 90° mirror which 
directed the light back to the pendulum. The beam then 
retraced its path to the autocollimator. A quarter-wave 
plate ensured that very little light bound for the detec- 
tor returned to the laser diode. The optical power in 
the light beam was ~ 15 juW. This optical system gave 
a detector signal that was nominally insensitive to the 
'swing', 'wobble' and 'guitar-string' motion of the tor- 
sion pendulum but had twice the optical gain compared 
to our previous [29] autocollimators. 

The position-sensitive detector (United Detector Tech- 
nology 1x3 mm model) was a 3-terminal, analog, 
resistive-division device. The location of the light spot 
(focused to a diameter of about 0.1 mm on the detec- 
tor) was inferred from the ratio of the difference A and 
sum E currents received at the two ends of the detector. 
The laser beam was chopped at 141.7 Hz, and Stanford 
Research Systems model SR830 lock-in amplifiers gener- 
ated the A and £ signals that were digitized and recorded 
by the data-taking computer. The time-constants of the 
lock-in amplifiers were normally set to nominal values of 
10 s. 



H. Vacuum Vessel 

The vacuum vessel was used in a previous short-range 
test of the equivalence principle [29]. A small ion pump 
kept the pressure inside the chamber near 10~ 6 Torr. The 
vacuum vessel was placed on a vibration-dampened kine- 
matic mount directly attached to an unused cyclotron 
magnet in our laboratory. The magnet was relatively 
insensitive to building vibrations and sway because it 
had an independent foundation and a mass of about 200 
tonnes. The vessel could be tilted with leveling screws to 
align the attractor assembly with the free-hanging pen- 
dulum as described in Sec. IV A 2 



a factor of ^80. The constant residual field inside the 
chamber was <5 mG. 



J. Environmental Monitors 

We continuously monitored the in-phase and quadra- 
ture outputs of the A and S lock-ins as well as several 
important environmental parameters: 

1. two Applied Geomechanics Instruments model 755 
electronic tilt sensors, attached to the top of the 
vacuum vessel, monitored apparatus tilt. 

2. vertical and horizontal Teledyne Geotech model S- 
500 seismometers, mounted on the vacuum vessel 
support frame, monitored seismic disturbances dur- 
ing Experiment II. 

3. eight Analog Devices model AD590 temperature 
sensors monitored the temperature variations of the 
apparatus and the air around it: 3 were mounted 
on the vacuum vessel, 4 measured air temperatures 
inside and outside the thermal shield, and 1 mea- 
sured the temperature of the radiator water bath. 
All sensors were encased in small copper blocks to 
increase the short-term stability of the readings. 

4. the ion-pump current served as a vacuum gauge. 



K. Calibration Turntable 

The gravitational calibration, described in Section 
V F 1 , employed an external turntable that rotated 2 
brass spheres uniformly around the vacuum vessel. The 
turntable was made from aluminum and driven by a mo- 
tor identical to that which drove the attractor and had a 
similar angle index marker. The turntable platter had 
6 evenly spaced holes sized to seat the two 5.08 cm- 
diameter brass spheres used in the gravitational calibra- 
tion. 



I. Mechanical Support and Shielding 

The entire apparatus and sensitive electronics were en- 
closed in a thermally insulating box made of 5-cm-thick 
foam. An external device (Neslab model RTE-221) circu- 
lated constant-temperature water through a radiator/fan 
assembly inside the thermal enclosure. The air tempera- 
ture inside the box varied by < 30 mK /day and variations 
at the 10o> signal frequency were less than about 1 mK. 

A single layer of 0.76-mm-thick mu- metal surrounded 
the lower half of the vacuum chamber. A second mu- 
metal plate was mounted between the motor and the vac- 
uum vessel. This shielding reduced the ambient field by 



L. Data Acquisition System 

We recorded data with a 80286 computer containing 
an Analogic LSDAS-16 interface card. The card had 16 
analog input channels with 16-bit analog-to-digital con- 
verters, one of which was multiplexed to record up to 
16 temperatures. Digital input registers counted the at- 
tractor drive pulses, giving us the instantaneous relative 
angle of the attractor, and the index pulse that gave us 
absolute angle information. Digital outputs controlled 
the attractor rotation direction and turned on and off 
the rotation. 

Data were taken at regular intervals determined by the 
computer's clock, typically between 5 and 15 s. The data 
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strobing interval was set to r /4n where r D was the nom- 
inal period of free torsional oscillations and n an integer 
that ranged between 8 and 17. 



IV. ALIGNMENT AND POSITIONING 

A. Positioning the Attractor Relative to the 
Electrostatic Shield 

We aligned the attractor with respect to the electro- 
static shield by temporarily replacing the Be/Cu foil and 
holder with a plate containing a hole through which we 
passed a depth micrometer probe to measure the distance 
to the attractor. This plate had the same thickness as one 
of the aluminum rings that held the electrostatic shield 
membrane as they were lapped simultaneously. By mea- 
suring the distance above each attractor leveling screw 
and making iterative adjustments, we made the attractor 
parallel to the lapped plate (and therefore to the electro- 
static shield) to within 0.1 mrad. The distance d am be- 
tween the membrane and the upper surface of the attrac- 
tor, 47±1 /im (54±3 fjxa) in Experiment I (II), was deter- 
mined by averaging the results of the depth-micrometer 
measurements, 48±3 /jm (57±3 /im), with values in- 
ferred from membrane-attractor capacitance measure- 
ments, 47±1 /im (51±1 /im). 



1. Leveling the Pendulum 

The pendulum ring was leveled (made perpendicular 
to the fiber) by using it as one plate of a differential ca- 
pacitor as shown in Fig. 6. The other two plates were 
semi-circular, horizontal copper pieces, separated by a 
thin gap, that were installed in place of the electrostatic 
shield. The pendulum was suspended ~ 1 mm above the 
plates and rotated uniformly using the fiber positioner 

9 stage. If the pendulum disk were perfectly horizontal, 
the differential capacitance between the two sides would 
not change as the pendulum was rotated. However, any 
pendulum tilt would cause the differential capacitance 
to vary at the rotation frequency. This Ilo signal was 
independent of any static tilt of the copper plates but 
was a strong signature of a tilted pendulum. Figure 7 
shows the differential capacitance signal for one revolu- 
tion of the stage before and after leveling. Initially, the 
differential capacitance was dominated by the Xlo signal. 
However, as the pendulum approached level, the Ilo sig- 
nal became so small that the modulation caused by the 

10 holes passing over gap between the plates became the 
dominant harmonic component (this 10w signal indicates 
that the gap between the semi-circular plates was not 
precisely aligned under the torsion fiber axis). We it- 
erated the leveling process until the lev differential ca- 
pacitance signal corresponded to a misalignment of only 
0.03±0.01 mrad (0.12±0.01 mrad). 




FIG. 6: Schematic diagram of the leveling measurement. Sup- 
pose the pendulum were tilted so that the "+" plate has a 
higher capacitance as shown in the left diagram. The right 
diagram shows the pendulum after the upper fiber attachment 
has been rotated by 180°; now the "— " plate has the higher 
capacitance. 




Angle (Degrees) 

FIG. 7: [color online] Differential capacitance for one revolu- 
tion of the pendulum before and after leveling. The tilt signal 
occurs at Ilo; the IQuj ripple results from the ten holes pass- 
ing over the gap between the plates of the capacitor. The ltd 
amplitudes correspond to pendulum tilts of 7.52±0.06 mrad 
and 0.03±0.01 mrad before and after leveling, respectively. 

A false leveling signal could occur if the fiber did not lie 
directly along the rotation axis of the stage. In this case 
the center of the pendulum ring would have moved in a 
circle with respect to the copper plates and the "run-out" 
would have modulated the differential capacitance at the 
rotation frequency. The fiber run-out was measured with 
a theodolite and, after adjustments of the upper fiber 
attachment, was reduced to 74±23 /im (24±8 /im) for 
Experiment I (II) . The effect of any stray signal induced 
by the residual fiber run-out (or stage and shaft motion) 
was checked as follows. First, we "parked" the pendulum 
on the capacitor plates and turned the 9 stage through 
~3 revolutions. Then, the pendulum was raised to the 
original height and left to unwind freely. With nothing 
but the pendulum's tilt to influence the differential ca- 
pacitance, any difference in the peak-to-peak amplitude 
between the freely unwinding case and the case when 
the stage was rotating must have been due to a false ef- 
fect. The observed difference introduced an additional 
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uncertainty of ±0.11 mrad (±0.01 mrad). The result- 
ing value for the pendulum tilt in Experiment I (II) was 
0.03±0.11 mrad (0.12±0.01 mrad). 



2. Aligning the Attractor 

Because our calculated Newtonian, Yukawa and power- 
law torques ignored most imperfections of the attractor 's 
rotary motion, we had to reduce these effects to negligible 
levels. We measured the attractor's wobble and vertical 
and horizontal runouts (see Figure 8) with high-quality 
mechanical indicators. By adding thin shims between the 
attractor and holder, we reduced the wobble amplitude to 
< 0.2 mrad (0.05 mrad) in Experiment I (II). We found no 
measurable vertical runout within the ±3 /jm measuring 
error. The horizontal runout was ±5 /im. 

We made the pendulum ring parallel to the conducting 
membrane using a capacitive technique. We tilted the en- 
tire apparatus and measured the pcndulum-to- membrane 
capacitance as a function of the tilt using a low-noise, 
high-sensitivity Stanford Research Systems model SR720 
LCR meter; because the capacitance is inversely propor- 
tional to the gap between flat plates, the minimum ca- 
pacitance occurs when the pendulum and membrane are 
parallel. The capacitance at fixed z-actuator setting was 
measured for several values of the tilt along two axes and 
fitted with a parabolic function of tilt angle. The results 
for one scan are shown in Figure 9. The tilt was then ad- 
justed to the capacitance minimum along each axis. The 
error on this procedure was ±18 /trad, determined by the 
accuracy of the tilt adjustment (±6 //rad) and the error 
on the fitted value of the minimum (±17 /irad). 



B. Centering the Pendulum on the Attractor 

The displacement between the pendulum and attractor 
was a key parameter in our experiments. The electronic 
actuators were sufficiently accurate (after accounting for 
the scale factors described in Section IV C) and repro- 
ducible to give the relative x, y, and z positions of the 
pendulum. It was more difficult to determine x a , y and 
z ; the x, y, and z values at which the pendulum would 
be resting on the attractor with its holes centered on the 
corresponding holes in the attractor. 



Wobble Rumble 




FIG. 8: [color online] Possible attractor rotation imperfec- 
tions. Vertical rumble was too small to measure with our 
instruments. 
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FIG. 9: Data and parabolic fit for the pendulum-to-membrane 
alignment along one tilt axis. The capacitance is minimum 
when the pendulum and membrane are parallel. The error on 
the leveling, ±18 /trad, was determined by the error on the 
fitted minimum position and the accuracy to which the tilt 
could be adjusted. 

1. Horizontal Centering 

Our horizontal alignment relied upon the gravitational 
torque on the pendulum being symmetric about x and 
y ■ We took data at points offset along the x and y axes 
and used these points as part of our complete data set, 
the uncertainties in the x and y values were taken to be 
±2 /tm, the specified reproducibility of the stages. The 
analysis left the alignment points, xq and yo, as uncon- 
strained free parameters. The fit was able to resolve the 
true alignment to ±10 /im. Because the torque varia- 
tion over an offset range of ±100 /im was much less than 
our statistical uncertainty, this alignment was more than 
adequate. 

2. Vertical Separation 

We extracted the absolute vertical location of the mem- 
brane from "z-scans" where we measured the pendulum- 
to-electrostatic-shield capacitance, C(z), as a function of 
z-actuator reading, z, using the Stanford Research In- 
struments SR720 LCR meter. We checked that the AC 
fields used to measure the capacitance did not apprecia- 
bly deform the thin shield membrane by taking data at 
drive voltages of 0.1 V and 1.0 V. 

We fitted the C{z) data with 

C(z)=MC th (z-z*) + C st , (12) 

where Af was a normalization constant that accounted 
for the absolute calibration of the capacitance meter, C s t 
was the stray capacitance of the suspension fiber, etc., 
and z* was the z-actuator reading corresponding to zero 
vertical separation between the pendulum and the elec- 
trostatic shield. The expected capacitance, Cth(C)i a t 
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TABLE IV: Comparison of the expected and fitted parame- 
ters of Equation 13 for the z-scan of Figure 10. 



Parameter 


Expected Value 


Fitted Value 


A/ 1 


1.00 ±0.01 


0.976 ± 0.005 


z* 


N/A 


1.067 ± 0.002 mm 




N/A 


0.78 ±0.11 pF 


s b 


~ 10 /xm 


5 ± 14 ^m 



a pendulum-shield separation £ was calculated with an 
accurate finite-element model of the pendulum and its 
electrostatic shield using ANSYS. 

The torsion fiber was a linear spring as well, so 
that seismic vibrations excited the pendulum's vertical 
"bounce" mode. Ultimately, this motion limited the min- 
imum achievable separation. We also needed to account 
for bounce when determining the membrane position. 
The bounce amplitude was measured by raising the pen- 
dulum until its mirror intersected only half of the auto- 
collimator light beam so that any vertical motion of the 
pendulum would alter the amount of light returned to 
the detector. By measuring the change in the amount 
of returned light for a known change in z we determined 
that AL/Az = 220 mV/mm, where AL and Az are the 
changes in returned light and height, respectively. We 
monitored the 11.5 Hz (9.7 Hz) bounce oscillation on a 
spectrum analyzer and found the average bounce ampli- 
tude to be 13 /im, with a maximum value of 36 /im over 
the course of one day. 

The act of moving the vertical stage during a z-scan 
excited the bounce mode. We modeled the bounce as 
a sinusoidal modulation of the pendulum-to-membrane 
separation so that Equation 12 was modified to 



C(z) = TV 



C th (z) + 



si d 2 c th (z) 



dz 1 



+ C st , (13) 



where <5^ is the mean square bounce amplitude. 

One z-scan typically contained about 25 different 
pendulum-shield separations between ~ 0.1 — 3 mm. The 
data points and fit for one z-scan are shown in Figure 10. 
The fitted values of the parameters in Equation 13 are 
compared to their best estimates in Table IV. The un- 
certainty in z» was typically ±3 /im. The z-actuator 
reading, z at which the pendulum would just touch the 
attractor was 



Zq — Z% t 



(14) 



where t is the membrane thickness and d am is defined in 
Section IV A. The separation between the bottom surface 
of the pendulum and the upper surface of the attractor 
is then 



(15) 




The uncertainty 5z = ±5 (im was found by combin- 
ing in quadrature the uncertainties in z* (±3 /im), t 
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FIG. 10: z-scan data and residuals from Experiment I. The 
fit indicated the pendulum would touch the membrane at a 
z-actuator reading of 1.067±0.002 mm. The lower panel dis- 
plays the fit residuals. The smaller set of error bars on the 
residuals represents the uncertainty in the capacitance values 
alone, while the larger ones include the effect of the ±0.5 jj,m 
uncertainty in the relative stage positions. 



(±3 fim) and d am (±2 /im); the uncertainty in the z 
readings was obtained by combining capacitance and z- 
micrometer data and varied between 3 (1.6) /xm and 6 
(4) /im in Experiment I (II). For the upper- attractor- 
disk-only runs in Experiment I the vertical position at 
vanishing pendulum-to-upper-disk separation, z' , was 
more uncertain, Sz' — ±0.025 mm, because c? am was 
measured with a less precise micrometer. 



C. Alignment and Calibration of the Translation 
Stages 

The alignment of the horizontal stages with respect to 
the plane of the electrostatic shield was checked by mea- 
suring the change in the pendulum-shield capacitance as 
the x and y-stages were moved (see Figure 11). We found 
a small misalignment (x-stage: 8 mrad, y-stage: 2 mrad) 
that was taken into account when computing the vertical 
separation values used in the data analysis. 

The absolute calibration of the 3 electronically- 
controlled stage actuators caused us much trouble. Ini- 
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FIG. 11: Dependence of the vertical separation on the hor- 
izontal position; the slope shows that the a;-stage was tilted 
by 8 mrad from the horizontal. 
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FIG. 12: Calibration of the electronic z-actuator with a pre- 
cision mechanical micrometer. The slope of the best fit line 
is 1.0230±0.0002, not the expected slope of 1.0000. The error 
bars are smaller than the size of the points. 



tially we assumed that the actuator readings were accu- 
rate to within their specifications. However, the Exper- 
iment I data did not agree with Newtonian predictions. 
Despite much effort, we could not find any systematic er- 
ror that could account for the results. Experiment II was 
designed solely to confirm this discrepancy, and showed 
a similar, but not identical, discrepancy. Finally, having 
exhausted all other possibilities, we decided to check if 
the digital actuators were correct by measuring the ac- 
tual distance moved by each digital actuator with preci- 
sion manual micrometers. The results are summarized in 
Figure 12 and Table V. We found highly linear relation- 
ships but the slopes were not unity! The position values 
used in our analysis were corrected with the appropriate 
calibration. 



TABLE V: Ratio of apparent displacement of the computer 
controlled actuators to the actual displacement determined 
using a precision manual counterpart. The ratios reveal that 
the stages did not move as far as expected. 



Axis 


Ratio 


X 


1.0239 ±0.0002 


V 


1.0248 ±0.0002 


z 


1.0230 ±0.0002 



V. TORQUE MEASUREMENTS 

A. Protocols 

Data runs, which lasted about one day, were charac- 
terized by definite values of the vertical and horizontal 
positions of the pendulum. The data-sampling interval 
and the attractor rotation rate were both synchronized 
to the pendulum's nominal free-oscillation period which 
was t = 399.7 s (467.0 s) in Experiment I (II). All sen- 
sors were recorded precisely 32 (68) times per free tor- 
sional oscillation period. The attractor was rotated with 
a period of r a tt = 17to (15to), so that the fundamental 
lOw signal repeated every 1.7 (1.5) torsion cycles; this 
placed the fundamental IQuj signal at a sufficiently high 
frequency to suppress the 1/f noise (see Appendix A), 
but kept the 20w and higher harmonics from being too 
strongly attenuated by pendulum inertia. We analyzed 
signals up to 30w (20a;); higher harmonics were attenu- 
ated below the autocollimator noise level. 



B. Extracting the harmonic twist amplitudes 

1. Converting the Autocollimator Output to a Twist Angle 

We converted the autocollimator difference A and sum 
S signals into twist angles using 



C-2 



(16) 



where ci and C2 are linear and quadratic calibration coef- 
ficients. The quadratic term arose from small differences 
in the electronic gains of the signals from the two ends of 
the photo-detector as well as from non-linearities in the 
autocollimator's position-sensitive detector. We deter- 
mined c 2 = —0.0416 by sweeping a light spot across the 
sensitive axis of the detector and minimized our sensitiv- 
ity to C2 by periodically adjusting the 6*-stage to compen- 
sate for fiber drift, keeping the mean position of light spot 
roughly constant. Our determination of ci is discussed 
in Sec. VF1. 
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2. Dynamic Calibration of the Autocollimator 

We repeatedly checked the linearity of the deflection 
scale and measured the free-oscillation and decay times of 
the torsion oscillator with dynamic calibrations. Energy 
was first taken out of the pendulum's torsional mode by 
manipulating the 9 stage until the pendulum had a very 
small twist amplitude, typically A ~ 5 /irad. Then data 
taking was started. At a time t Q after beginning the run 
the 9 stage was abruptly rotated by a known amount, 
typically A$ = 349 ^trad (0.02 degrees). The calibration 
coefficients were checked by fitting the pendulum's ob- 
served response 9{t) (which was an implicit function of 
ci and C2) to the expected change 

9^(t) = a + /3i+L4sin(fi* + 7 )+ (17) 

e(t - t„)A0(i - cosn(t - * ))]e-(*-*°)/ Td , 

where a and /3 account for fiber "drift" , A and 7 ac- 
count for the initial oscillation of the pendulum, fl is the 
free-oscillation frequency, is a unit step function, and 
Td is the decay time of the pendulum amplitude. Fig- 
ure 13 shows the data from one such calibration. The 
dynamic calibration uncertainty was dominated by the 
±8.5 fiY'&d (2.5%) uncertainty in the (9-stage position 
(half the minimum angular step size and resolution of 
the encoder) . Dynamic calibrations determined the nom- 
inal free-oscillation period, r Q = 2tt/CI, and exponential 
amplitude-damping time, Td (needed for the inertial cor- 
rections described in Section VE) and provided an ex- 
cellent test of the linearity of the gross angular deflection 
scale. However, because the amplitude of the dynamic 
calibration data was necessarily much larger than the 
that of our ISL data (~ 5^irad) we did not use these data 
for the final calibration of the deflection (and therefore 
the torque) scale, but instead applied a known gravita- 
tional torque as described in Sec. VF. 



3. Filtering out the Free Oscillations 

The top panel of Figure 14 shows typical raw A au- 
tocollimator data. The free torsional oscillation was the 
dominant component of the pendulum's twist (we de- 
liberately kept this amplitude large to facilitate precise 
measurements of 51). Because this motion provided no 
information about an ISL violation, we applied a digital 
filter to suppress the free oscillations before proceeding 
with the analysis. This "notch" filter averaged two data 
points separated by one half of a torsional period, and 
assigned the average value to their midpoint. 



(4>) 



1 r 
2 



('-t) + 8 ('+t) 



(18) 



where t is the time when the attractor was at angle 0, 
t = 2ir/D, is the torsional period, and 9 is the interpo- 
lated value of the twist. The middle panel of Figure 14 
shows the filtered results. 
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FIG. 13: Dynamic calibration data. At t=580 s, the upper ro- 
tation stage changed the equilibrium position by 0.349 mrad. 
The dashed line depicts the equilibrium position as a function 
of time. The fit (not shown) passes so close to the data points 
that it cannot be distinguished from them in this plot. 



4- Harmonic Analysis of the Twist Data 

Each data run was subdivided into separate "cuts" 
containing exactly two oscillations of the lOw signal (72° 
of attractor rotation). For each cut, 9f was fitted as a 
function of attractor angle <f> with 



'th(^) = ^2[bn sm(n<p) + c n cos(n^)] + 



m=0 



d P 



(19) 

All points in a given cut were treated equally and ar- 
bitrarily assigned an uncertainty ±1 ADU. The d m co- 
efficients of the Legendre polynomials P m accounted for 
fiber drift. The first sum in Eq. 19 ran over n =1, 3, 
10, 20, and 30 (the n =30 term was not needed in Ex- 
periment II). The floating b n and c„ coefficients for the 
n = 10, 20, and 30 terms yielded the "gravitational" sig- 
nals. The n = 3 coefficients accounted for the three-fold 
symmetry of the attractor drive system while the n = 1 
coefficients allowed for any effect periodic at the attrac- 
tor rotation frequency. (When fitting the individual cuts, 
the small n = 1 and n — 3 terms were fixed to values ob- 
tained from fitting data covering one complete attractor 
revolution.) The smooth curve in the middle panel of 
Figure 14 shows the fit and the bottom panel shows the 
fit residuals. For each cut we computed a "chi-square" 
of the A^j data points 



x 2 



=X> f (&)-M&)) a , 



(20) 



that reflected the goodness of the fit. 
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FIG. 14: Top panel: raw autocollimator data for one attrac- 
tor revolution (five cuts) at a vertical separation of 234 /im 
in Experiment I. Middle panel: the same data after it has 
been filtered and fitted with Eq. 19. The ten oscillations re- 
flect the tenfold rotational symmetry of the instrument. The 
asymmetric shapes are due to phase shifts in the twist-to- 
torque conversion; the IOoj frequency is less than Q while the 
20a) and 30u) frequencies are greater than fl. Bottom panel: 
the fit residuals. 
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FIG. 15: Upper panel: pendulum twist signal; lower panel: 
vacuum chamber pressure. 



of a total of 1599 (676) cuts In Experiment I (II). Tests 
described below showed that the noise in the remaining 
cuts was normally distributed. 



1. Pressure Spikes 

The first rejection criterion was an abrupt disturbance 
of a sensor. The only observable irregularities were spikes 
in the vacuum chamber pressure which often visibly ex- 
cited the pendulum (see Figure 15). We rejected all cuts 
with pressure spikes greater than ten standard deviations 
above the mean pressure. This criterion rejected 108 (5) 
cuts in Experiment I (II). 



C. Rejecting Bad Data 

We excluded cuts from our final analysis that were pol- 
luted by external disturbances with nongaussian distri- 
butions. These "bad cuts" were identified by a spike in 
the vacuum chamber pressure, an anomalously large \ 2 , 
or an abrupt change in the amplitude or phase of the free 
torsional oscillation. These criteria rejected 228 (27) out 



2. Bad x 2 of the Harmonic Fit 

An unusually large value of x 2 indicated that the pen- 
dulum had been disturbed in some way. The x 2 distri- 
butions and rejection levels for Experiments I and II are 
shown in Figure 16. This criterion rejected 111 (14) of 
the cuts passing the first criterion in Experiment I (II). 
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FIG. 16: Distributions of x 2 values for all 1599 (676) cuts 
in Experiment I (II). The dashed vertical lines show the cut 
rejection levels. The rightmost bins contain all cuts with x 2 's 
extending past the ends of the plots. 



Jumps in the Free Torsional Oscillation Amplitude and 
Phase 



The final rejection criterion was an abrupt change in 
the amplitude or phase of the free torsional oscillation. 
We determined the two quadrature amplitudes of the 
free oscillation in each cut and found the differences be- 
tween adjacent cuts Aa S m and Aa cos - The difference 



J = ^/Ac 



Aaj? os is the "jump" in the free oscilla- 



tion amplitude. If Aa S i n and Aa cos are distributed about 
zero with a common spread, and there are M different, 
normally-distributed sources producing jumps, each with 
its own average amplitude <r, the jump distribution has 
the form 



dN(J) 



M N- 

2—1 L 



- j2 /^dJ , (21) 



where Ni are the numbers of events of each type. Fig- 
ure 17 shows the jump data from both Experiments. 
Equation 21, with M = 2, and N and a% as free parame- 
ters, provides an excellent fit to the observed jump distri- 
butions, giving x 2 ' l v °f 0.98 (0.68) in Experiment I (II); 
on the other hand, fits with M = 1 were poor. The 
two stochastic sources may, for example, be vertical and 
horizontal seismic vibrations to which the apparatus re- 
sponds quite differently. We eliminated cuts with jumps 
that were not distributed normally by setting the rejec- 
tion level at the point where the fits began to deviate 
significantly from the data: 2.4 /irad (1.2 /xrad) in Exper- 
iments I (II). The rejection level was lower in Experiment 
II because the magnetic damper, which did not function 
properly in Experiment I, was fixed before beginning Ex- 
periment II. This criterion rejected 9 (8) of the surviving 
cuts in Experiment I (II). 



D. Combining Good Data 

All N c cuts within a given run that survived the rejec- 
tion criteria were treated as equivalent so that the mean 



FIG. 17: Distribution of jumps in the 2 experiments. The 
curves, based on Eq. 21 with M — 2 provide good fits to the 
data to the left of the vertical dashed lines, indicating that 
these jumps are normally distributed with two independent 
sources with <ti = 0.32 ± 0.02 ^rad (0.20 ± 0.02 /xrad) and 
(72 = 0.78 ±0.05 ^rad (0.41 ±0.03 ^irad) in Experiment I (II). 
The vertical dashed lines are the rejection levels; larger jumps 
are due to spurious events such as small earthquakes, pressure 
spikes, or spontaneous unwinding of the torsion fiber. Note 
that in thermal equilibrium the mean free oscillation ampli- 
tude is 9 rulB = ^/ks T/k = 1.1 (ir&d, where ks is Boltz- 
mann's constant. 



and its standard deviation of the Eq. 19 coefficients from 
an individual run were 



1 Nc 



(i) 



\ N C (N C 



1 



/](K(i) ~ b r , 
t=l 



(22) 



(23) 



similar equations apply to the c„ coefficients. Figure 18 
shows that the b n and c„ coefficients of the individual 
accepted cuts were essentially normally distributed. 



E. Converting Twists to Torques 

The applied torque on the pendulum, N, at frequency 
uj s was inferred from the twist angle 9(lo s ) using 



(24) 



where I is the rotational inertia of the pendulum, N and 
9 are complex numbers whose real and imaginary parts 
represent the components varying as cosw s t and sinw s t, 
respectively; / = /i/t c /f is a complex quantity that ac- 
counts for the effects of pendulum inertia, electronic time 
constants and digital filtering — its imaginary part corre- 
sponds to a phase shift. We extracted the torques using 
Eq. 24 rather than N — k6 (where k is the torsional 
spring constant of the fiber) because electrostatic cou- 
plings (see Section VI A 4) can cause the effective k to 
vary slightly with s, whereas / is truly constant. 
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The differential equation describing the twist of a pen- 
dulum subjected to a driving torque N and to anelastic 
damping (which acts like an imaginary component of the 
torsion constant [36]), is 



16 + k 1 



N 



(25) 



where Q = is the quality factor of the torsion oscil- 

lator. A periodic driving torque, N = Noe luJt , produces 
a steady-state twist 



0(t) = 



1 



N 



i - (uj/n) 2 + i/Q J m 2 



(26) 



where il = -J k/I. Equation 26 shows that /; = 
A\ exp where 

CI 2 



A i = V[l-(w/fi) 2 ] 2 + (l/Q) 2 ; ^i^tan- 1 



Q(n 2 -cu 2 ) ' 

(27) 

The torsion filter described in Section V B 3 has a re- 
sponse /f = Af expi$f, where 



. 



(28) 



In computing the effect on 9 of the digitally synthesized 2- 
pole filters of the A and £ lock-in amplifiers, we could ne- 
glect the time constant of the £ amplifier because its sig- 
nal was nearly constant in time. Then / tc = A tc exp z<I>tc, 
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FIG. 18: Normalized residuals n = (6j — bi) / (\f~NiO-i) for 
the accepted data, where en is given by Eq. 23. Top panel: 
all 5574 IOcj, 20u> and 30oj b n and c n coefficients of the 2-disk 
attractor data from Experiment I. Bottom panel: all 2596 Wcu 
and 20lu coefficients from Experiment II. The smooth curves 
are Gaussian fits that give a = 0.90±0.01 and a = 0.91±0.01 
respectively, indicating that our assigned errors (Eq. 23) were 
slightly conservative. 
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FIG. 19: Response of the lock-in amplifier after a steady input 
was abruptly removed at time to ■ The fitted curve (largely in- 
distinguishable from the points in this plot) assumes an ideal 
2-pole filter response with time constant t tc , where to and t tc 
are free parameters. 



where 



A tc = (1 + ^ 2 4) ; $tc = 2 tan 1 (cut tc ) 



(29) 



The lock-in amplifiers had nominal time constants of 
10 s; we found actual value of ttc by using a function 
generator to drive the amplifier's analog and reference 
inputs with sine-wave and square-wave signals, respec- 
tively. The time dependence of the output after the in- 
put was abruptly disconnected, shown in Fig. 19, revealed 
that t tc = 10.082 ±0.011 s . 



F. Calibrating the Torque Scale 

1. Absolute gravitational calibration of Experiment II 

The torque scale was calibrated in Experiment II by 
applying a known gravitational torque on the pendu- 
lum. Two 7.938±0.003 mm-diameter aluminum spheres 
with average masses (corrected for buoyancy in air) 
of M s = 0.7348 ± 0.0001 g were placed in diametri- 
cally opposed holes on the pendulum ring, producing 
known qi m gravitational inner multipole moments, where 
I and m define the spherical multipole order. At the 
same time, two brass spheres with average masses (in 
air) of 515.1809 ± 0.0016 g[30] were placed in diametri- 
cally opposite holes in the calibration turntable platter 
13.983 ± 0.004 cm from the rotation axis. These spheres 
produced known Qi m outer multipole moments so that 
as the calibration turntable rotated at a frequency u> c , it 
produced a torque on the pendulum given by[29] 



N, 



cal 



8ttG 



^ 21 

1=2 



I 

tE 

m=2 



lm{mqi m Qi m e imu " t ) ■ (30) 
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FIG. 20: Data and fit for three revolutions of the calibration 
turntable. The 2u c signal is dominant, as expected. 



FIG. 21: Absolute calibration data as a function of vertical 
separation. The curve is a polynomial fit whose curvature 
agrees well with the multipole prediction. 



The sum was terminated at I — 8 because higher-order 
torques were negligible. We considered only the m = 2 
component of the torque, which was calculated from 
the masses and locations of the spheres to be iV™i =2 = 
40.042 ± 0.039 fN-m. Figure 20 shows typical data. Be- 
cause the 13.98 cm length scale of our calibration was 
close to the 16.76 cm scale used in measuring[31] G, any 
ISL violation (see Fig. 1) could not have significantly ef- 
fected the accuracy of our calibration. 

We fitted data with the calibration turntable displaced 
from its nominal center position by a series of horizontal 
and vertical points and corrected for the small measured 
backgrounds from mass asymmetries in the pendulum 
body and the calibration turntable; the maximum twist 
was 6»™f 2 = 12.6330 ± 0.0048 ^rad, Figure 21 shows the 
observed signal as a function of vertical position and the 
best fit based on the multipole prediction (the effects of 
displacements on the qi m and Qi m moments were calcu- 
lated using multipole translation theorems [29]). 

The torque corresponding to a twist signal t9 s was com- 
puted as 



tcrminc 



N(uj s 



i n 2 /M) 

Ica,l ^al /(2w C) ^ca] 



(31) 



cal 



where / = 175 ± 1 g cm 2 is the normal pendulum's rota- 
tional inertia, I ca i = I + / sp h where I sp h = 10.544 ±0.004 
g cm 2 is the rotational inertia of the small calibration 
spheres used on the pendulum. The uncertainties in /, 
I sp h, N™= 2 , ^cai and 6>™" 2 produce a torque scale- factor 
error of 0.12%. We folded in an additional 0.5% error to 
account for autocollimator non-linearities, measured in 
data taken subsequent to Experiment II, to give an over- 
all torque scale factor C\\ — 1.000±0.005; this was incor- 
porated into our analysis as described in Section VIII A 
below. 

Note that the left hand side of Eq. 31 is independent of 
the autocollimator coefficient c\ defined in Eq. 16. How- 
ever we can use the gravitational calibration data to de- 



ci = 1.2758 ± 0.0073 mrad 



(32) 



2. Relative calibrations in Experiment II 

We made relative gravitational calibrations of each 
Experiment II data run to test for and correct any z- 
dependence of the torque scale. The calibration turntable 
with the two brass spheres was rotated continuously at 
frequency lo c . The qn moments of the pendulum frame 
and Q14 moments of the brass spheres on the calibration 
turntable interacted to produce a 4w c torque. The cali- 
bration turntable was raised and lowered along with the 
pendulum to ensure the calibration torque was indepen- 
dent of z. Figure 22 shows the observed 4uj twist as a 
function of vertical position. The small z-dependence of 
the signal could have been caused by irregularities in the 
reflectivity or curvature of the pendulum mirror. The 
torques for each point were corrected for this variation of 
the calibration. 



3. Calibration of Experiment I 

Because Experiments I and II used the same autocol- 
limator and suspension fiber, we could infer the calibra- 
tion for Experiment I from the absolute calibration of 
Experiment II. We adopted the c\ value determined by 
the Experiment II gravitational calibration (Eq. 32) but 
increased the uncertainty to ±1.5 % to account for the 
uncertainty in / and autocollimator non-linearities. Ex- 
periment I torques were then computed using Eq. 24; we 
adopted a scale-factor uncertainty C\ = 1.000 ± 0.015. 
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FIG. 22: Calibration torque signal (relative units) as a func- 
tion of the pendulum's vertical position in Experiment II. 



4- Checking the calculated rotational inertias 

The torques measured using Eq. 24 were proportional 
to the pendulum rotational inertias which were calculated 
from detailed models of the pendulums. We verified the 
calculation of / for Experiment II using the known rota- 
tional inertia of the spheres / sp h and the measured free- 
oscillation frequencies from the gravitational calibration 
data 

Q = 13.4584 ± 0.00027 mrad/s (33) 
ftcai = 13.0709 ± 0.00173 mrad/s , (34) 

where fi was extracted from data with s > 370 /im. We 
accounted for the effect on the torsional constant k from 
fiber stretch due to the extra weight of the spheres 



Ak _ ^Ad 
K d 



AL 
~L 



(35) 



where d w 20 /im and L w 800 mm are the diameter and 
length of the fiber. We determined AL from the pendu- 
lums's observed bounce frequency and total mass M 
so that 



AL 
~L 



2M s g 



{2wu h ) 2 ML 



= -0.00012 



(36) 



The tabulated Poisson's ratio for tungsten gives Ad/d = 
-0.28AL/L, so that Ak/k = -2.12AL/L, or 



1.00026 



Therefore 



hi + hph 



K 
K cal 



n 2 



(37) 



(38) 



which can be solved to give In = 175.3±0.8 g-cm 2 , which 
agrees well with the value 175 ± 1 g-cm 2 computed using 
our detailed model of the pendulum. 



TABLE VI: Harmonic torque components from the on-center 
data in Experiment I. All torques have units of fN-m. The 
errors do not include the overall 1.5% scale-factor uncertainty. 



s (mm) 




Nio 








N20 






JV30 






0.216 


5 


327 ± 





051 


2 


358 ± 


040 


2 


866 ± 





031 


0.234 


5 


283 ± 





022 


2 


287 ± 


015 


2 


794 ± 





024 


0.382 


4 


471 ± 





032 


1 


998 ± 


026 


2 


344 ± 





036 


0.577 


3 


729 ± 





044 


1 


671 ± 


031 


1 


835 ± 





026 


1.067 


2 


051 ± 





041 


1 


046 ± 


031 


1 


054 ± 





033 


1.556 





974 ± 





052 





642 ± 


043 





608 ± 





046 


2.045 





140 ± 





030 





406 ± 


021 





387 ± 





035 


3.022 


-0 


752 ± 





040 





228 ± 


023 





139 ± 





037 


3.999 


— o 


970 ± 


o 


034 


o 


059 ± 


031 





045 ± 





033 


4.977 


-1 


031 ± 





026 





066 ±0 


028 


-0 


027 ± 





038 


6.443 


-0 


909 ± 





042 





063 ±0 


028 





008 ± 





031 


8.132" 


5 


005 ± 





038 





035 ±0 


033 





030 ± 





075 


8.134 at 


5 


008 ± 





065 





077 ±0 


396 





056 ± 





228 


8.134 a 


4 


993 ± 





046 





067 ±0 


028 





036 ± 





042 


8.798 a 


3 


846 ± 





066 





075 ±0 


048 





044 ± 





076 


8.817 a 


3 


823 ± 





057 





035 ±0 


043 





044 ± 





037 


10.773 a 


1 


780 ± 





045 





058 ±0 


027 





039 ± 





025 


10.773 a 


1 


830 ± 





025 





050 ±0 


022 





063 ± 





033 



"Upper-disk-only runs 

^Run taken at twice the attractor rotation frequency of the other 
runs 



We checked the calculated / for Experiment I in a sim- 
ilar fashion, using the measured free-oscillation frequen- 
cies of the Experiment I and II pendulums and correcting 
for the effect on k of the differing weights of the two pen- 
dulums using equations similar to Eqs. 36 and 37. We 
found h = 128.5 ± 0.7 g-cm 2 which agreed well with the 
detailed model calculation which gave 128 ± 1 g-cm 2 . 



G. Torques from Experiment I 

Experiment I contained 27 runs, each of which con- 
sisted of at least 36 cuts. The data were collected over a 
period of ~ 45 days. Several runs were taken with the up- 
per attractor disk alone at vertical separations s > 8 mm. 

The full range of vertical separations was 216 /mi to 
10.77 mm, including the upper-disk-only runs. The mea- 
sured harmonic components of the torque as a function 
of s are shown in Table VI. N w changes sign at a separa- 
tion of ~ 2 mm where the torque from the lower attractor 
plate overcomes that from the upper plate. Additional 
runs were taken at varying horizontal offsets to determine 
the unknown alignment point (x ,y ). 



H. Torques from Experiment II 

Experiment II contained 17 runs, each of which con- 
sisted of least 19 cuts. The entire data set was collected 
over a period of ~ 30 days. The full range of vertical sepa- 
rations was 137 /im to 7.90 mm. The measured harmonic 
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TABLE VII: Harmonic torque components for the on-center 
data runs of Experiment II in units of fN-m. The errors do 
not include the 0.5% scale-factor uncertainty. 



s (mm) 


Nio 


N20 


0.137 


31 287 ± 081 


10.024 ± 0.042 


1 ^7 

U. J-O / 


OL.A'-±\J ZIZ U.VJ i I 




0.162 


30.750 ± 0.094 


9.842 ± 0.034 


0.191 


30.079 ± 0.072 


9.625 ± 0.029 


0.372 


26.382 ± 0.047 


8.368 ±0.023 


0.763 


19.671 ±0.052 


6.169 ±0.028 


1.545 


10.030 ± 0.037 


3.349 ±0.021 


3.011 


1.008 ±0.023 


1.057 ±0.018 


3.700 


-0.992 ±0.016 


0.655 ±0.018 


5.459 


-2.828 ±0.015 


0.183 ±0.014 


7.903 


-2.492 ±0.018 


0.056 ±0.012 



components of the torque as a function of separation are 
shown in Table VII. In this measurement the Nio phase 
reversal occured at a separation s ~ 3.3 mm. 

Six of the 17 runs were taken for horizontal centering 
and were not used to constrain ISL violation because the 
relative calibration setup was not employed. These 6 runs 
were used only to determine the central values and errors 
of x and y . 



VI. SYSTEMATIC ERRORS 

A. False Effects 

False effect are any effects, other than the gravitational 
interaction between the pendulum ring and the attractor 
disks, that produced a real or apparent coherent torque 
on the pendulum at our signal frequency. Whenever pos- 
sible, we tested the sensitivity of our apparatus to indi- 
vidual false effects by exaggerating the effects and report 
the systematic effect on the fundamental 10ui signal. Un- 
less otherwise noted, the effect on the higher harmonics 
was negligible. We took explicit account of those sys- 
tematic effects that were not negligible compared to the 
typical 8(4) nrad statistical error per Experiment I(II) 
run. 



1. Thermal 

We measured the temperature-twist coupling by sub- 
jecting the apparatus to a Ilo temperature modulation 
and measuring the pendulum's response (the large ther- 
mal mass of the apparatus prevented us from obtain- 
ing an adequate temperature variation at higher frequen- 
cies). Three sensors attached to the vacuum vessel quan- 
tified the temperature change. The radiator water bath 
had a luo temperature amplitude of 910±1 mK, while the 
average of the three sensors mounted on the vacuum ves- 
sel recorded a 123±0.5 mK \uo variation. The induced 



twist amplitude was 1.17±0.07 /irad, corresponding to a 
9.51±0.57 /irad/K temperature feed-through. The vari- 
ations in the temperature at the fundamental lOw fre- 
quency were < 0.1 mK during normal data taking so that 
the systematic error associated with temperature modu- 
lations was <1 nrad. 

A false signal from thermal effects could also have oc- 
curred because heat generated by the drive motor made 
the attractor slightly higher warmer than its surround- 
ings. When the motor drive was energized, its tempera- 
ture increased by 1 K. Consider Experiment I. The holes 
in the upper attractor plate were essentially blocked off 
by the lower plate, so that information about the attrac- 
tor motion could be transmitted thermally only to the 
extent that the top surface of the upper attractor plate 
was at a different (lower) temperature than the top sur- 
face of the lower plate. This temperature difference AT U \ 
could have produced rotating "hot spots" on the elec- 
trostatic shield that might have transmitted information 
about the attractor rotation to the pendulum. 

Several time scales are relevant for assessing the size 
of such effects. The time scale for the foil temperature, 
T to relax by radiation is 



7"rad 



c pP t 

8<ieT 3 



55 s 



(39) 



where p ~ 8900 kg/m 3 is the density of copper, C p and p 
are the heat capacity at constant pressure and density of 
copper, <7 is the Stefan-Boltzmann constant, and e ~ 0.1 
and t — 20 ^m are the emissivity and thickness of the 
foil, respectively. The time scale for the foil temperature 
to relax by conduction through the copper is 



C pP L 2 
4k 



0.2 s 



(40) 



where k is the thermal conductivity of copper and L, 
taken to be 1 cm above, is the distance between 2 points 
on the foil. The time scale for the upper attractor disk 
temperature to relax by residual gas conduction to ob- 
jects at a temperature T, separated by less than a mean- 
free-path, is 



C p pw 



7psJk B /(mT) 



3 x 10 5 s , 



(41) 



where ks is Boltzmann's constant, m is the mass of a 
gas molecule (assumed to be H 2 ), p <~ 10~ 3 Pa is the 
pressure, and w = 2 mm is the thickness of the upper 
attractor plate. 

We estimated A7^i by noting that the lapped plates 
are in good thermal contact and that the time constant 
for the top plate to cool via radiation (~ 10 5 s) is much 
longer than the time constant for thermal conduction be- 
tween the plates so that we can assume the bottom sur- 
face of the upper plate is at the same temperature as 
the top surface of the lower plate. The energy lost by 
the upper plate through radiation to the foil is supplied 
by thermal conduction through the upper plate, giving a 
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temperature difference between the top and bottom sur- 
faces of the upper attractor plate, which is essentially 
AT ul , 



This leads to a torque 



AT ul = 



4eVuT 3 AT uf 
k 



where ATf is temperature difference between the upper 
attractor plate and the conducting foil, and e' = e/2 ~ 
0.05 was obtained by assuming that the upper plate and 
foil have the equal emissivities. Taking an upper limit 
AT uf ~ 1 K, wc find that AT ul ~ 1.5 x 10~ 6 K. 

The temperature variations on the foil (AT in Eqs. 46, 
47 and 48) due to ATi can be estimated as follows. The 
rotation period of the attractor, ~ 7000 s, is much longer 
than the thermal equilibration time of the foil so that we 
can consider the case of a stationary source. If we model 
the thermal variations in the azimuthal coordinate of the 
foil as a one-dimensional thermal conductor subject to 
a spatially-sinusoidal heat source, then the steady-state 
heat equation becomes: 



d 2 T _ S cos (2irx/L) 
dx 2 k 



(43) 



where S is the rate of heat input per unit volume. The 
solution to this equation gives us AT = SL 2 /(4ir 2 k) — 
e'aL 2 T 3 AT ul /(2Tr 2 kt) where t = 20 fim is the thickness 
of the foil. Taking L ~ 2 cm as the distance between the 
holes in the source plate, we find 



AT ~ 6 x 10 



-10 



K 



(44) 



Finally, we estimate the effect on the pendulum of hot 
spots on the foil as follows. Assume that moving spots, 
which are hotter by an amount AT, are approximated as 
squares of side length 2a where a is the radius of a hole. 
The outgassing rate q(T) of the foil depends strongly on 
the temperature T, 



g(T) = Qexp(-e/T) , 



(45) 



where is an activation temperature that has a char- 
acteristic value in the neighborhood of 10 4 K. Assuming 
that the pressure, p, is dominated by the outgassing rate, 
Eq. 45 implies that a local temperature difference AT 
will produce a local pressure difference 



Ap 

P 



6AT 
~T 2 ~ 



(46) 



The molecules travel ballistically because the mean free 
path is long compared to the relevant dimensions. The 
maximum force on a pendulum hole occurs when an at- 
tractor hole is displaced horizontally by a. Molecules 
emitted by a hot area 2a 2 hitting an area ah, where h is 
the thickness of a pendulum hole, apply an extra hori- 
zontal force on a single hole 



AT 



ahp 9AT 



(47) 



.e 



N ~ 10 3 — AT fN-m , 



(48) 



(42) 

where we assumed p <~ 10 3 Pa. Equations 48 and 44 



predict a thermally induced torque is 4 x 10 fN-m or 
a spurious signal of ~ 0.01 nrad. 

We checked for a temperature gradient effect by de- 
liberately heating the attractor with resistors attached 
to the attractor drive motor, raising the temperature of 
the motor mount by 4.8 °C — a five-fold increase over the 
temperature rise associated with turning on the motor. 
The 10w twist signal in this configuration was consistent 
at the 1.5cr level with the predicted negligible size of this 
effect. 

Thermal expansion from the heating associated with 
turning on the motor could have affected the attractor- 
to-membrane distance. The attractor-membrane capac- 
itance was unchanged within ±0.01 pF out of 1.33 pF 
when the motor was turned on, corresponding to a dis- 
tance change of < 1 /im for Experiment I. We therefore 
assumed that systematic errors associated with the ther- 
mal gradient produced by the motor were negligible. 



2. Tilt 

Tilt can produce two separate false effects. The first 
arises if the apparatus tilt is modulated at the signal fre- 
quency Tilting the fiber attachment is known to twist 
the pendulum, so that any tilt at the signal frequency 
would produce a false signal. We tilted the vacuum ves- 
sel by 1.02 ±0.01 mrad and found that the pendulum 
equilibrium position changed by 16 ± 2 /xrad. We ana- 
lyzed the tilts recorded by the AGI sensors in each nor- 
mal run and an upper limit of 4.0 nrad on signal periodic 
at 10w; this corresponded to a tilt systematic at lOw of 

< 0.06 nrad. 

The second tilt effect arises if there is a static tilt be- 
tween the pendulum and attractor. Our calculations of 
the expected Newtonian and Yukawa torques assumed 
that the plates were perfectly parallel. In normal data 
taking, the tilt between the pendulum and attractor was 

< 0.2 mrad (see Section IV A 1). The effect of ignor- 
ing this tilt was calculated as follows. Let the calculated 
lOw gravitational signal of an untilted pendulum at a 
pendulum-to-attractor separation of s be N\q(s). A tilt 
of size tp « 1 would change the effective separations 
between the centers of masses of the holes. This effect 
causes the torque on a tilted pendulum 



10 

£ 



N w (s + r p V> sin(j7r/5)) 
10 



(49) 



to depend quadratically on ip (the sum runs over the 10 
pendulum holes and r p is the distance of the holes from 
the pendulum's symmetry axis). The maximum expected 
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difference for any Experiment I (II) run was 

IJVin — JV tnt I — 1 o w m-4 /t a w in-4 



ln | - 1.2 x 1CT 4 (1.6 x 10" 4 ) fN-m . (50) 



Tilt also causes the direction of the symmetry axes of 
pendulum holes to deviate from the assumed vertical. 
But this second-order effect is smaller than that given in 
Eq. 50 by an amount ~ (a p /r p ) 2 and therefore can be 
neglected. 

Before beginning Experiment I we checked the static 
tilt effect by deliberately tilting the apparatus by 1 mrad 
and measuring the change in iVio- We attempted to rc- 
center the pendulum after the tilt but uncertainties in 
the centering prevented us from being able to verify the 
very small effect predicted by Eq. 50. We adopted the 
value in Eq. 50 for the systematic error due to tilt. 



3. Magnetic Background 

We measured the magnetic moment of the pendulum 
by applying a time-varying external magnetic field. A 
coil placed outside the thermal shielding produced a hor- 
izontal 6.3 mG field at the pendulum location. We 
modulated this field at lOw and saw induced twists of 
335 ± 79 nrad (252 ± 25 nrad) . Under normal operating 
conditions, the amplitude of lOw magnetic field varia- 
tions at the site of the pendulum, measured by a flux-gate 
magnetometer, was < 0.6 [iG. Hence, the error associated 
with a magnetic coupling was ±0.03 nrad (±0.02 nrad), 
much smaller than our statistical errors. 



AU « ae^Vp^A/ (2dp m ) where A is the total area of the 
10 pendulum holes and d pm ~ 60/im is the minimum 
pendulum-to-membrane separation. The moving bumps 
therefore exert a torque on the pendulum 



N = AU/A<j) < 



which is negligible. 



[ Km " 


2 




1 mV_ 




1 mV 



2 x 10~ 15 fN-m 



(52) 



5. Gravitational Background 

Our calculations of Newtonian and Yukawa torques as- 
sumed that the torque was due entirely to the holes in 
the pendulum ring and attractor disk. However, resid- 
ual gravitational coupling between the pendulum frame 
and the attractor drive mechanism could have produced 
small spurious signals at lOw, 20w, etc. The apparatus 
was designed to suppress such couplings — the pendulum 
frame had two-fold (four-fold) rotational symmetry and 
the attractor drive was 3-fold symmetric — so we did not 
expect significant torques from these components at any 
harmonic of interest. We tested for residual torques at 
the signal frequencies by taking data with the attractor 
disks removed from their holder and found no resolved 
effect. We then replaced the attractor in its original po- 
sition and removed the pendulum ring from its frame and 
again found no signal. Because we did not expect or ob- 
serve any gravitational background signals, we assigned 
a negligible systematic error for this effect. 



4- Electrical Background 



6. Electrostatic Effects on k 



Suppose the attractor, membrane and pendulum were 
ideal conductors but at slightly different potentials with 
voltage differences V am and V mp . The attractive pressure 
between two flat conducting plates with a potential dif- 
ference V and separated by a distance d is e V 2 /(2d 2 ). 
This pressure is essentially absent over the holes, so the 
membrane feels a missing pressure above each attrac- 
tor hole P am = e V 2 m /(2dl m ), where d am w .05 mm 
(see Sec. IV A). This will excite the m = 10, n = 1 
drumhead mode of the membrane (the lowest mode that 
could produce a lOto electrostatic coupling of the attrac- 
tor to the pendulum). The deflection in this mode is 
S(r, 4>) — aJio(Zcior) cos(10<^)), where Jio is a Bessel func- 
tion of the first kind and 



a = 



C P a m 

4tV/i 2 „,i 



1 mV 



6 x 10~ 16 mm 



(51) 



C w 0.07 is the projection of the attractor hole pat- 
tern on the shape of the membrane mode, pt is the areal 
density of the membrane, and its resonant frequency is 
/io,i = 6.02/o,i ~ 5.1 kHz (see Sec. HID). These rotat- 
ing "bumps" will in turn change the electrostatic bind- 
ing energy of pendulum-membrane system by an amount 



Direct electrostatic interactions between the pendulum 
and attractor were essentially eliminated by the conduct- 
ing membrane. However, real conductors differ from ideal 
ones and electrostatic potential wells on the surface of 
the membrane could interact with patch charges on the 
pendulum to change the effective k of the torsion fiber, 
making it a function of z — z* . Figure 23 shows the tor- 
sion period, 2tt/Q as a function of separation for Experi- 
ments I and II. A small systematic variation in the period 
occurred at very close separations. We sidestepped this 
problem by measuring the free oscillation frequency f2 for 
each run, and inferring the torques from the measured 
twist using Eq. 24 which does not explicitly depend on 
k. Note that any effect that changed the oscillator prop- 
erties was eliminated by this procedure. 



7. False Effect Error Budget 

Table VIII summarizes the systematic twist error pro- 
duced by false effects. The sum of these effects is far 
smaller than our statistical noise, which was ~ 8 (4) nrad 
per data point. 
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FIG. 23: The torsion period as a function of the pendulum- 
to-membrane separation in Experiments I (upper panel) and 
II (lower panel). The horizontal lines show the mean periods 
for separations > 100 /im). 



TABLE VIII: Contributions to the systematic twist error from 
dominant false effects. 



Effect 


lOw 


Twist Error (nrad) 


Temperature Modulation 




<1 


Apparatus Tilt 




<0.06 


Relative Tilt 




<0.05 


Magnetic 




<0.03 



B. Errors in the Twist-to- Torque Conversion 

We checked the validity of the frequency-dependent 
twist-to-torque conversion described in Section VE by 
comparing data taken at different attractor rotation rates 
u). In Experiment I, where the normal rotation period 
was T a tt = 17t , we also took data with attractor pe- 
riods of 7t and 43t . The results are summarized in 
Table IX. As we found no measurable discrepancy and 
expected none, we assigned a negligible error to the con- 
version (except for that introduced by measurement un- 
certainties in Q, t d , uj and t tc ). 



TABLE IX: Comparison of the torques measured at two dif- 
ferent attractor rotation speeds in Experiment I. 



Attractor Period 


Measured Torques [fN-ml 




lOw 


20a; 30uj 


7t 


3.204 ±0.040 


1.564 ±0.116 1.630 ±0.144 


43to 


3.227 ±0.037 


1.590 ±0.019 1.702 ±0.052 


Signal change 


0.023 ±0.054 


0.026 ±0.118 0.072 ±0.153 



C. Uncertainties in the Instrumental Parameters 

Uncertainties in important instrumental parameters 
(such as the pendulum-to-attractor separation, torque 
scale calibration, and the missing masses of the holes) 
affect the predicted Newtonian and Yukawa torques. We 
measured these parameters and used their values and er- 
rors as constrained parameters when fitting the torque 
data, in effect letting gravity give us an independent mea- 
sure of these experimental parameters. The parameters 
were C (the overall torque-scale factor), M p (the total 
missing mass of the pendulum holes), M" p (the mass of 
the upper attractor holes), M x ° (the mass of the lower 
attractor holes) , z D (the vertical position reading at van- 
ishing pcndulum-to- attractor separation), g (the gap be- 
tween the upper and lower attractor disks), and ip (the 
angular misalignment of the two attractor disks). In ad- 
dition for Experiment I we needed z' Q (the z value for the 
upper-disk only runs) and for Experiment II we needed 
M" p ' out (the mass of the out-of-phase holes in the upper 
attractor) and ip' (the misalignment angle between the 
in-phase and out-of-phase holes in the upper attractor 
plate). 



1. Missing Mass of the Holes 

The missing masses of the holes were key factors in 
our experiment. In Experiment I, we inferred the missing 
masses from the densities of the materials and the hole 
dimensions. The densities were determined after we had 
taken all the data for that experiment. We measured the 
masses of the pendulum ring and the attractor disks, and 
then drilled 10 extra holes in the pendulum ring and 8 
extra holes in each of the attractor disks. These masses, 
together with the measured volumes of the normal and 
extra holes gave us the missing masses of the normal 
holes. 

A more direct technique was used to determine the 
missing masses of the holes in Experiment II. Each disc 
was weighed before the active holes were drilled (Mi) 
and then again immediately after drilling the holes (M2). 
Since the upper plate had 2 sets of holes, the mass was 
measured after drilling each set. The attractor disks were 
then lapped to the final thickness and weighed again 
(A/3). (The pendulum was not lapped so that, for it, 
M 3 = A/2.) Then the gold layer was applied and the 
piece weighed again (A/4). The resulting hole masses 
were 

Mi 

M = (Mr - M 2 ) -f + / Au (M 4 - M 3 ), (53) 

where the computed factor /a u accounted for the gold 
missing from the flat surfaces of the holes as well as that 
present on the round edges of the holes. 
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constraint on g to minimize the effect of correlations be- 
tween g and other fitting parameters that would increase 
the uncertainties in the extracted physics parameters a 
and A. 

In Experiment I, the upper and lower disks had thick- 
nesses of 1.846 ±0.005 mm and 7.825 ±0.005 mm, re- 
spectively. When the disks were stacked the combined 
thickness was 9.681 ± 0.007 mm, yielding g = (0.010 ± 
0.010) mm. We used the same technique in Experi- 
ment II, but took many more thickness measurements 
and used a more precise micrometer to obtain g — 
(0.0014 ± 0.0010) mm. 



FIG. 24: Radial positions of the individual inner holes of the 
Experiment II upper attractor disk. 
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FIG. 25: Residuals of the angle between the individual inner 
holes of the Experiment II upper attractor disk. 



4- Angular Misalignment of the Attractor Disks 

The angular alignment of the upper and lower attractor 
disks was maintained by three pins that passed through 
holes in the upper and lower disks. We determined the 
angular misalignment, ip between the two disks by ex- 
amining the disk assembly with a measuring microscope. 
We found that if = 0.00±0.01 degrees (0.13±0.01 de- 
grees) for Experiment 1 (II). 



VII. PREDICTING THE EXPECTED 
NEWTONIAN, YUKAWA, POWER-LAW AND 
MASSIVE PSEUDOSCALAR TORQUES 



Conventions 



2. Hole Placement 

The calculated torques were based on the average val- 
ues of the distance of the holes from the rotation axis 
and the angle between them. It was important to know 
whether there were any systematic deviations from the 
average values (for example, the holes could have been 
on an ellipse instead of a circle). Figures 24 and 25 show 
the radii and angles of the inner holes of the upper at- 
tractor disk in Experiment II, respectively. There are no 
systematic variations or outliers. The holes in the other 
disks (and pendulum rings) used in the two experiments 
displayed similar behavior. 



We assume that the pendulum ring and attractor disks 
are perpendicular to the vertical torsion fiber, so that a 
torque about the fiber axis can only be caused by a hori- 
zontal force F t , perpendicular to the fiber and tangential 
to the pendulum ring. The horizontal force between any 
two perfectly circular holes is directed along the horizon- 
tal projection, t, of the line between their centers, so that 
the torque exerted by the pendulum on the fiber is 

N = r{i) x F t (i,j) + £ f(i) x F t (i, k) , (54) 

i,j i,k 



3. Gap between the Attractor Disks 

Any gap g between the attractor disks (for example, 
because they were not perfectly flat) would affect the 
cancelation of the 10w torques. Although g was a pa- 
rameter in fitting our data, we needed an independent 



where i runs over the 10 pendulum holes, j and k run over 
the holes in upper and lower attractor disks, respectively, 
and F t has Newtonian, Yukawa, power-law and massive 
pseudoscalar components 

F t = F t G + a(\)F?(X) + l3 k F t p (k) + 7 (A)F t M (A) . (55) 
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where a, (3k and 7 arc dimcnsionlcss parameters. The 
forces F t G , , F t p and F t M were computed as 



F t G (t,s) 



F?(t,.s,\) 



F t p (t,s,k) 
F t M (t,s,X) 



Gpip 2 
Gpip 2 
Gpip 2 
Gpip 2 



VtJV-2 



dV^S 

at 

dV 1 dV 2 ^- 



dVidV 2 



vJv. 



dt 

d_ 

dt 



1 

R 

exp(— R/X) 



(56) 



v 1 Jv 2 



dVidV 2 ^- 
at 



R 

r \ k j_ 

R 

r \2 K^R/X) 



A 



where R = \R 2 — R\\ is the distance between infinitesi- 
mal volumes in the two cylinders, ro = 1 mm, and the 
integrals extend over the volumes V\ and V 2 of the two 
cylinders. 




attractor 
rotation 
axis 



FIG. 26: Geometry used to calculate the torque on the pen- 
dulum from the rotating attractor. The attractor hole is on 
the left. 



B. Calculating the Harmonic Torque Database 

For the geometry shown in Figure 26, the torque about 
the torsion hber axis from a single pendulum/attractor 
hole pair is 

N^a, S, <f> 8 ,8) = r p F t (t,s) (57) 
/ r a sin (f> a + S sin <f> s \ 
X \ t ) ' 

where F t (t, s) is the total horizontal force on the hole at 
a pendulum-to-attractor vertical separation s, and 



t = y (r a sin 4> a + S sin (f>s ) 2 + (r a cos 4>a + 6 cos <j>s - r p ) 2 . 

The procedures for efficiently calculating F t G , F t y (A), 
F t p (fc) and F t M (A) are described in Appendix B. 

The torque on the pendulum from one complete set of 
attractor holes contains 100 terms, 

JV(0 B ,M«,s)= (58) 
9 9 

^^^1(^ + ^0+^/5, 6, (j> s -kTr/5,s) , 

j=0 fc=0 

where the actual attractor angle is (f> a — (f> + (j), with <p 
being the attractor angle reading and <p is a (f> value at 
which the torque vanishes. We computed ./V for cj)$ = 
and <ps — 7r/2; because the results were identical to 
within 10~ 18 N-m, it was sufficient to fix 4>$ = 0. Fig- 
ure 27 shows the calculated torques from the upper and 
lower disks of Experiment I for attractor angles between 
— 7r/10 and n/10 and a vertical separation of s = 0.2 mm. 

The tenfold rotational symmetry of the attractor was 
exploited to compute a database of the harmonic torque 
amplitudes at frequencies nui, 

20 r /w ~ 

N n {6,s) = — N((f> 1 0,5,0,s)sm{n<j))d<j) . (59) 
Jo 
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FIG. 27: Calculated Newtonian torque on the Experiment I 
pendulum at s = 0.2 mm as a function of attractor rotation 
angle. The dash-dotted and dotted curves are the torques 
from the upper and lower attractor disks, respectively. The 
solid curve is the total torque magnified by a factor of ten. 



for 20 different S's between and 2 mm and 100 values of 
s. Numerical evaluation of this integral is exact as long 
as a sufficient number of equally spaced points are used. 



C. Generating Interpolated Torque Functions 

We generated two-dimensional analytic harmonic 
torque amplitudes for each set of attractor holes at each 
signal frequency of interest by cubic-spline interpolation 
of the harmonic torque database. Each potential required 
six interpolating functions: one for each of the 2 (3) sets 
of holes and the 3 (2) frequencies in Experiment I (II). 
Figure 28 shows the interpolated lOw torque from both 
attractor disks of Experiment I for 5 = as functions of 
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FIG. 28: [color online] Predicted Newtonian IOcj torque on 
the Experiment I pendulum as a function of s. The solid 
curve shows the total torque multiplied by a factor of 10. 
The dash-dotted and dotted curves show the torques from 
the upper and lower disks, respectively. 

vertical separation s. 

D. Precision Requirements 

We wanted the overall absolute precision of the torque 
calculations, SN, to be much better than the statisti- 
cal uncertainty of the experimental values, which was as 
good as SN w 1 x 10~ 17 N-m (see Table VII). Because of 
the large cancelation between the torques from the up- 
per and lower attractor disks, SN « y/2 6N&, where SN^ 
is the precision of the calculated torques on the pendu- 
lum from an individual upper or lower attractor disc. We 
therefore require 

5N d « 7 x 1(T 18 Nm . (60) 

The torque Nd = 10rF t d where r = .03 m and the force 
on a single pendulum hole, F^, consists of a sum over 10 
attractor holes. Although Fig. 29 shows that the force on 
a given pendulum hole is dominated by the forces from 
the two nearest attractor holes and that these two forces 
tend to cancel, SF? m VlO SF t . We satisfy Eq. 60 if 

5FP » -^L- « 7 x 10- 18 N . (61) 

We adopted the requirement SF^ < 1 x 10~ 19 N. Two 
sources of error had to be considered. The desired preci- 
sion determined the stopping point of the numerical inte- 
gration procedure, as well as the spacing of the database 
points that were used to interpolate the force at arbi- 
trary s and t. Figure 30 shows the calculated Newtonian 
force and the errors associated with the database and 
its interpolation. The contributions to the error of the 
Newtonian torque calculations are listed in Table X. 
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FIG. 29: The solid line shows the horizontal Newtonian force, 
Ft 3 , °n a single pendulum hole in Experiment I (scaled to a 
maximum value of unity) as a function of horizontal hole sep- 
aration t. One attractor hole has been placed where the force 
between cylinders is maximum (t ~ 5.4 mm). The horizon- 
tal distances of the other 9 attractor holes are shown. The 
vertical lines show the contributions of each upper attractor 
hole to the force vector. The horizontal dashed line is the net 
force from all 10 attractor holes; it is largely determined by 
the two attractor holes closest to the pendulum hole. 



TABLE X: Error budget for the Newtonian torque calcula- 
tions. 



Source 


Error 




Force Calculation 


< 1 X 10" 2U 


N-m 


Force Interpolation 


< 1 x 10" 20 


N-m 


Torque Interpolation 


< 2 x 10" 20 


N-m 


Total 


~ 2.5 x 10~ 20 


N-m 



Our calculations of the Yukawa, power-law and mas- 
sive pseudoscalar forces could be less precise because, 
in essence, these forces were compared to the Newto- 
nian residuals of our data that were always within about 
2<T of zero. Because the fractional uncertainties in the 
best- fit values of a, 0% and 7 from uncertainties in the 
measured torques were greater than 0.4, the systematic 
uncertainties from the corresponding calculations of the 
forces between the two closest holes (see Fig. 29) could 
be neglected as long as 

5 Ft 

"TT « 0.4 . (62) 
-ft 

We computed the non-Newtonian forces with an abso- 
lute accuracy good enough to satisfy Eq. 62 in all cases. 
For example, the fractional accuracy of the Yukawa force 
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both experiments were fitted by minimizing 



FIG. 30: A typical Newtonian force calculation. The heavy 
line shows Ft from the upper attractor disk of Experiment I 
at a separation of s = 250 /im. The light line shows the in- 
tegration error, determined by doubling the number of steps 
in evaluating each of the two dimensions of the force integral. 
The points show the error associated with the spline interpo- 
lation which we found by comparing the values interpolated 
halfway between the data-base points to a direct numerical 
evaluation of the force integrals. 



calculation was ~ 4.5 x 10 2 at A = 10 /xm, 
A = 100 /xm, and 2 x 10~ 4 at A = 1 mm. 7 



10" d at 



VIII. RESULTS 



A. Fitting Functions 



We constrained ISL violations by comparing the mea- 
sured complex nuj torques, N n (Q) ± 5N n (Q), at pendu- 
lum positions indicated symbolically as Q — (xj,yj, Zj), 
to calculations of the expected Newtonian and possi- 
ble Yukawa, power-law, or massive pseudoscalar torques. 
The calculated torques were explicit functions of Q as 
well as of P key instrumental parameters (the missing 
masses M p , M^ p , Mj? of the holes in the pendulum and 
the lower and upper attractors, the alignment parame- 
ters x , j/o, z , 0o, and the angular offset <p and gap g 
between the two attractor disks) denoted symbolically by 
the array r) m , and the ISL- violating parameters a, A and 
Pk denoted symbolically as <;. The complete data sets for 
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(63) 



where j refers to the particular run, n to the signal fre- 
quency (lOw, 20w, etc.), and m to the individual exper- 
imental parameters (P max of which were independently 
constrained) with rf^ and <5?/°* p being their indepen- 
dently measured values. For Experiment II we included 
an uncertainty in the phase of measured torques, nd(j} , 
where 8<fi = ±0.046 degrees was deduced from the devi- 
ation of the index angle readings. We accounted for the 
uncertainty in the independent variable s by setting 

A7V„ = yf(6N n )* + (6s dN n /dsf . (64) 

Because x , y and 4> were not independently measured 
they were not included in the last sum in Eq. 63. 

The two-disk attractor data of Experiment I were fitted 
with 

N n (0 , 77, = e m0 ° M p [M^N up (n, S 3 ,z 3 —z ,<;) (65) 
+ e mtp Ml°N lo {n, S^Zj-Zc + hi+g, <r) . 

where Af up and A^ (which were calculated assuming unit 
mass for the holes) have the forms 

N(n,5,z,<;) = N G (n,S,z) + aN Y (n,5,z,X) (66) 
+ l3 k N p {n, S, z, k) + 1 N M {n 1 8, z, A) , 



and 5j = \J (xj — xq) 2 + (yj — yo) 2 . The upper-disk- 
only runs of Experiment I used only the first term on 
the right-hand side of Eq. 65. 

The parameter M p , which multiplies all torques, is de- 
generate with the experimental torque scale factor C . We 
therefore assumed C = 1 exactly and combined the scale 
factor uncertainty with the measurement uncertainty in 
M p when specifying SM p . The Experiment II data were 
fitted with a function similar to Eq. 65 but with addi- 
tional terms that accounted for the out-of-phase holes in 
the upper attractor disk. 



B. Comparison of Data with Newtonian 
Predictions 

We first compared our data with the Newtonian pre- 
diction by fixing a, the /3fc's and 7 in Eq. 66 to zero. 



1. Experiment I 



The fit for Newtonian gravity alone was excellent with 



122, v = 156 and Q( X 2 , v) 



0.97, where 



27 



TABLE XI: Comparison of the measured and fitted experi- 
mental parameters for Experiment I. The uncertainty in M p 
reflects the uncertainty in the overall torque-scale factor as 
well as in the hole masses. Primed quantities refer to the 
upper-disk-only runs. 



Parameter 


Measured Value 


Fitted Value 


M p 


3.972 ± 0.060 g 


4.096 ± 0.014 g 


M a up 


11.770 ± 0.004 g 


11.7707 ± 0.0039 g 




88 fififi + 019 q- 


88 665 ± 019 p- 


XO 


N/A 


-0.413 ± 0.007 mm 


<;<> 


N/A 


-0.977 ± 0.012 mm 


Zo 


0.000 ± 0.005 mm 


0.001 ±0.002 mm 


4>o 


N/A 


0.007 ±0.007 deg 


x'o 


N/A 


-1.114 ±0.100 mm 


y'o 


N/A 


-0.846 ±0.100 mm 


i 

z a 


0.000 ±0.025 mm 


0.025 ±0.013 mm 




N/A 


-0.003 ± 0.023 deg 


g 


0.010 ±0.010 mm 


0.023 ±0.001 mm 


v> 


0.00 ±0.01 deg 


0.0095 ± 0.0006 deg 
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FIG. 31: [color online] Centering data and Newtonian fits as a 
function of the magnitude, S, of the offset. This plot includes 
points displaced in both x and y. Solid, dashed and dotted 
lines denote iVio, 2 x N20 and 1.5 x A30, respectively. 



Q(x 2 > v ) is the probability that a \ 2 as poor as that ob- 
served could arise by chance. The best fit values of the 
experimental parameters were all within about 2 stan- 
dard deviations of their measured values. Table XI dis- 
plays the measured and fitted values of the parameters. 
Figure 31 shows the Newtonian fit to the centering data 
at an average vertical separation of 234 jxva.. The top 
panel of Fig. 32 shows the data and Newtonian fit for 
all of the on-center runs; the bottom panel shows the 
residuals. 
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FIG. 32: [color online] Upper panel: torques and Newtonian 
fit for on-center data from Experiment I. Errors are smaller 
than the point size. Square, round, and diamond-shaped 
points show iVio, A20, and A30 torques. Triangular points 
are N10 of the upper-disk-only runs; these show the degree 
to which our 2-disk attractor attenuates the Newtonian sig- 
nal. Lower panel: residuals; the solid and dot-dashed curves 
show the expected N\q residuals from an a = 1, A = 250 /im 
Yukawa and from a /3s = 0.005 power-law interaction, re- 
spectively. The corresponding Yukawa curves in Figure 3 of 
Ref. [7] and Figure 2 of Ref. [6] were incorrectly designated 
as a — 3. 



2. Experiment II 

The Newtonian fit for the Experiment II data was 
good, with x 2 = 40.5 for v = 41, with Q(x 2 , v) = 0.40. 
Again all fitted values of the parameters were within 
about 2 standard deviations of their measured values. 
Table XII displays the measured and fitted values of these 
parameters. Figure 33 shows the data and Newtonian fit 
and residuals for all 11 on-center runs. 



3. Combined Results 

We fitted both data sets simultaneously with a func- 
tion consisting of the sum of the Experiment I and 
Experiment II functions. The combined fit assuming 
Newtonian gravity had \ 2 = 163.4 with v = 197 and 
Q{x 2 : v ) = 0.95. The values of the fitted and measured 
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FIG. 33: [color online] Top: Experiment II on-center data 
and Newtonian fit; solid square and round points show iVio 
and N20 torques, respectively. Bottom: residuals of the New- 
tonian fit. The uncertainties are larger at small separations 
(large signals) because the error contributions from the mea- 
sured free torsion period and in dN/ds have a bigger effects 
at small s. 



C. Constraints on Short-Range Yukawa 
Interactions 

We generated constraints on Yukawa ISL violations us- 
ing a fitting function that included the Newtonian and 
Yukawa terms. We analyzed each experiment individ- 
ually as well as the combined data set. As the results 
for Experiment I and Experiment II were consistent, we 
quote only the constraints extracted from the combined 
data sets. Data were fitted for 24 fixed values of A be- 
tween 0.020 mm and 10.0 mm, and for each case we ex- 
tracted the central value ot\ and its l-a uncertainty, Sax ■ 
We then used these values to establish 95% confidence 
intervals («2<r limits) on the allowed values of +a, —a, 
and |ev|. The constraints on \a\ were found by integrat- 
ing a normal distribution, N(a\, 5a\) with mean a\ and 
width Sa\, and finding the value of \ct\\ that satisfied 

/-|oaI poo 
N(a\,5a\)da\ + / N(a\, 5a\)da\ = 0.05. 

(67) 

Table XIII shows the 95% confidence intervals on pos- 
itive and negative a as well as the limit on |cv| for the 24 
different values of A, and Figure 34 displays the \a\ con- 
straint. It is easy to understand the general form of our 
constraints. In the limit where A is small compared to 
d, the characteristic dimension of the holes (in our case 
their thickness), the holes essentially behave like infinite 
parallel plates for which the Yukawa potential energy per 
unit area is 

— = 2ira\ 3 G Pl p 2 e~ s/x , (68) 



parameters for the combined data set were essentially 
unchanged. 



where s is the pendulum-to-attractor separation The 
torque is proportional to the derivative dV/ d(f> so that for 
a given torque sensitivity and smallest separation s m i n , 
we expect the constraints to have the asymptotic form 



TABLE XII: Comparison of the measured and fitted experi- 
mental parameters for Experiment II. The uncertainty in M p 
reflects the uncertainty in the overall torque-scale factor as 
well as in the hole masses, tp' is the angular offset between 
the in-phase and out-of-phase holes in the upper attractor. 



Parameter 


Measured Value 


Fitted Value 




2.6618 ± 0.0140 g 


2.6423 ± 0.0057 g 


M up,in 


8.6228 ± 0.0004 g 


8.6227 ± 0.0004 g 


M up,out 


19.2592 ± 0.0005 g 


19.2592 ± 0.0005 g 


M lo,in 


125.0984 ± 0.0006 g 


125.0984 ± 0.0006 g 




N/A 


0.803 ± 0.003 mm 


yo 


N/A 


-0.189 ± 0.003 mm 


Zo 


0.000 ± 0.005 mm 


-0.001 ± 0.003 mm 




N/A 


-0.001 ± 0.012 deg 


g 


0.0014 ±0.0010 mm 


0.0008 ± 0.0010 mm 


<p 


0.13 ±0.01 deg 


0.118 ±0.006 deg 


f 


0.00 ± 0.01 deg 


0.003 ± 0.009 deg 



loga(A) oc -3 log A ± s min /A ; (69) 
and an asymptotic slope on a log(cv) vs. log(A) plot of 
d log a s min 



d\og A 



A 



(70) 



which is satisfied by our constraints shown in Fig. 34. 
On the other hand, for A > > d, the cancelation from the 
lower plate reduces the signal from the Yukawa force and 
the constraints eventually weaken with increasing A. 

Our constraints improved upon the previous results 
shown in Fig. 1 by a factor of up to ~ 10 , and on more 
recent results [32, 33] shown in Fig. 34 by up to a factor 
of almost 10 2 . Our experiments are the only tests to date 
that reach gravitational sensitivity for length scales less 
than 500 fim. In particular, Yukawa interactions with 
\a\ > 1 are excluded at 95% confidence for A > 197 /im. 
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X [m] 

FIG. 34: [color online] Yukawa constraints from our com- 
bined data set as well as from other work[9-12, 32, 33]. The 
area above the heavy curves is excluded at the 95% con- 
fidence level. Predicted ISL violating effects from "extra 
dimensions" [1] , from dilaton[I5], moduli[16] and radion[I8] 
exchange, and from a conjectured solution[4] of the cosmolog- 
ical constant (vacuum energy) problem are shown as fainter 
lines. 



TABLE XIV: 68% confidence constraints on power-law po- 
tentials of the form given in Eq. 2 from this work and from 
previous work tabulated in Ref. [34]. 



k 


|/3fc|(this work) 


|/3fe| (previous work) 


2 


3.6 x 10" 3 


1.3 x 10 _3 [9] 


3 


2.8 x 10" 3 


1.3 x 10~ 2 [9] 




2.9 x 10" 3 


1.3 x 10 _1 [10] 


r, 
o 


2.3 x 10" 3 


2.1 x 10 _1 [10] 


TABLE XV 


: 95% confidence upper bounds on |7(A)| where 


A = h/(mc) 


and m is the pseudoscalar 


mass. 


A [mm] 


mc 2 [meV] 


It(a)| 


0.02 


9.85 


6.2 x 10 s 


0.05 


3.94 


2.6 x 10 3 


0.10 


1.97 


2.8 x 10 1 


0.20 


0.985 


1.24 


0.50 


0.394 


8.7 x 10~ 2 


1.0 


0.197 


2.5 x 10" 2 


2.0 


0.0985 


1.1 x 10" 2 


3.0 


0.0657 


8.0 x 10" 3 


5.0 


0.0394 


6.4 x 10" 3 


10.0 


0.00197 


5.5 x 10" 3 



E. Constraints on Couplings of Massive 
Pseudoscalars 



TABLE XIII: 95% confidence level constraints on Yukawa in- 
teractions from the combined data set. 



A (mm) a \a\ 

0.010 (-4.7 ±6.4) x 10 y 1.0 x 10 lu 

0.025 (-7.4 ± 10.4) x 10 4 1.6 x 10 5 

0.050 (-2.2 ± 7.9) x 10 2 8.8 x 10 2 

0.10 (0.2 ± 1.8) x 10 1 1.8 x 10 1 

0.25 (0.9 ±4.0) x 10 _1 4.3 x 10 _1 

0.50 (1.0 ±4.5) x 10" 2 4.8 x 10~ 2 

1.00 (0.1 ± 1.1) x 10" 2 1.1 x 10 -2 

1.50 (-1.8 ±7.2) x 10" 3 7.9 x 10~ 3 

2.50 (-5.1 ±5.9) x 10" 3 1.0 x 10~ 2 

5.00 (-7.3 ±6.7) x 10 -3 1.3 x 10~ 2 

10.0 (-0.7 ± 19) x 10" 3 1.8 x 10~ 2 



D. Constraints on Power-Law Interactions 



We constrained power-law violations of the ISL by fit- 
ting our combined data set with a function that contained 
the Newtonian term and a single power-law term. This 
procedure was carried out for power-law potentials with 
k = 2, 3, 4, and 5. The results are listed in Table XIV 
together with constraints from previous ISL tests given 
in Ref. [34]. 



Second-order exchange of massive pseudoscalars was 
constrained by fitting the combined data set with a func- 
tion containing the only the Newtonian and massive 
pseudoscalar terms. The results are listed in Table XV. 



IX. SOME IMPLICATIONS OF THE RESULTS 

A. Extra-Dimension Scenarios 

The most basic constraint from this work is an upper 
limit on the maximum size of an extra dimension. If we 
assume that one extra dimension is much larger than all 
the others, and compactification on a torus, ISL violating 
effects will first show up as a Yukawa interaction with 
a = 8/3 (see Eq. 3). Our 95% confidence level constraints 
require that 

i?» < 160 fim . (71) 

The scenario of two equal-size large new dimensions, 
whose radius i?* = 0.38 mm (see Eq. 4) is appropriate 
to give M* = 1 TcV, is ruled out by our results. Ac- 
cording to Eq. 4, our 95% confidence level constraint, 
X(a = 16/3) < 130 /im implies 

M* > 1.7 TeV/c 2 . (72) 
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B. New Boson- Exchange Forces 

1. Yukawa Interactions: Radion and Dilaton Constraints 

We provide significant constraints of the radion- 
mediated force predicted by Eqs. 6 and 7. For n = 1 
we find A(a = 1/3) < 270 while for n = 6 we have 
A(a = 3/4) < 215 //m which gives the following 95% 
confidence level constraints on the unification mass 



M*(n = 1) > 3.0 TeV/c 
M*(n = 6) > 4.6 TeV/c 2 



(73) 



Note that the radion constraint given in [33] is not valid 
as it was based on an earlier, erroneous prediction for a. 

Our 95% confidence level constraint on the dilaton 
mass, based on the couplings predicted in Ref. [15], is 



m > 1.0 meV/c 



2. Yukawa Interaction from Axion Exchange 



(74) 




g n 2 /Withc) xlO 7 

FIG. 35: [color online] Constraints on the couplings of mass- 
less pseudoscalar particles to protons and neutrons. The 
shaded region is excluded at 68% confidence. The down- 
sloping line is from this ISL test; the two curved lines are 
extract from a previous Equivalence-Principle test [29]. 



The predicted Yukawa interaction from axion exchange 
(Eq. 8) is too small to have been detected in this work. 
Even if the axion exists and has a mass in the upper por- 
tion of the allowed range (greater than 0.1 meV/c 2 ), its 
Yukawa interaction is unlikely to be detected in the near 
future. The axion's monopole-dipole and dipole-dipole 
interactions (discussed in Ref. [19]) seem more promising 
avenues for detecting axions in this mass range. 



3. Power-law Interactions 

The 75 couplings of massless pseudoscalar bosons to 
neutrons and protons, given by Eq. 11, are constrained 
by our limits on 03 from Table XIV. The contribution 
from coupling to electrons can be ignored because of the 
very small upper limit on such couplings deduced from an 
electron-spin-dependent experiment (see Ref. [34]). Con- 
sider the dimcnsionless parameter 



~z~ 




~Z~ 


4 


~N~ 




~N~ 


5n 




l 




2 (M 2 


. M_ 


l 


. M . 


2 (M 2 



(75) 
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where [Z/fj] and [N/fj] are defined in Ref. [29] and the 
subscripts 1 and 2 refer to the pendulum and attractor 
materials, aluminum and copper, respectively, with 



= 0.48181 



Al 



= 0.45636 



Cu 



Al 



= 0.51887 , (76) 
0.54475 . (77) 



Cu 



Our experimental limit on T, 
cC 

T = /3 3 — 64tt 3 u 4 (1 mm) 2 = 2.56 x 10" 10 /3 3 , (78) 

h A 

leads to the constraints on and shown in Fig. 35. 
Figure 35 also shows constraints we extracted from a pre- 
vious Equivalence Principle test [29] using Eq. 32 of that 
work[35] and assuming Aa — 2.8 x 10~ 13 cm/s 2 . 



4- Couplings of Massive Pseudoscalars 

Our ISL constraints on g^ and g 2 , shown in Fig. 35, 
apply to massless pseudoscalars. The corresponding con- 
straints on the couplings of massive pseudoscalars (see 
Eq. 10) are obtained using 

T(A) = 5.12 x 10- 10 7 (A) , (79) 

and can be obtained using Table XV or Fig. 36. 

C. Vacuum-Energy Scenarios 

The constraints on these initial experiments were not 
sensitive enough to check the suggestion [4, 5] that gravity 
essentially "shuts off" (i.e. a w —1) at Aw 100 /im. 

D. Future Prospects 

We are now running a second-generation instrument 
based on the same general principles but with a larger 
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FIG. 36: Scaling factor by which the 68% confidence upper 
bounds on g% and g\ from this experiment (shown in Fig. 35) 
should be multiplied to account for a non-zero pseudoscalar 
mass m. 

number of smaller-sized holes in the pendulum and at- 
tractor. This device has several other technical improve- 
ments that increase the signal and reduce the noise. We 
expect this new instrument will have the sensitivity to 
check the vacuum-energy scenario. The results will be 
reported in a later publication. 



to viscous drag, eddy currents etc., while the loss angle <f> 
accounts for internal friction[36] of the suspension fiber. 

In the familiar case of pure velocity damping (b > 0, 
4> = 0), the spectral density of torque noise is indepen- 
dent of frequency 

«H)=4fc B T^ , (A2) 



where f2 = y 1 njl is the free resonance frequency and Q = 
Ifl/b is the quality factor of the oscillator. This torque 
noise implies a spectral density in the angular deflection 
noise of 

(M-)> - Q j {n 2_ UJ 2 r + {nuj/Q) 2 ■ (A3) 

In this case there is no advantage, from the point of view 
of thermal noise, for placing the signal at any particular 
frequency. 

On the other hand, for pure internal damping (6 = 0, 
4> > 0) the spectral density of thermal noise has a 1/f 
character 

m 2 

(A&M) = 4fc B T— , (A4) 
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APPENDIX A: THERMAL NOISE 

Thermal noise in any oscillator sets a fundamental 
limit on the achievable statistical error on its amplitude. 
A single-mode torsion oscillator subject to both velocity 
and internal damping obeys the equation 

N = 10 + b$ + k(1 + i<f>)6 , (Al) 

where N is the applied torque, I the rotational iner- 
tia, 9 the angular deflection of the oscillator, and n the 
torsional spring constant of the suspension fiber. The 
velocity-damping coefficient b accounts for any losses due 



where now (5 = 1/0. The corresponding spectral density 
of thermal noise in the twist is 



{ th ^' QloI (f> 2 - c^ 2 ) 2 + (f! 2 /Q) 2 ' ( ' 

so that it is advantageous to boost the signal frequency 
above f2 until the noise contribution from the angular 
deflection readout system, {0 2 o }, becomes comparable to 

K). 

These points are illustrated in Fig. 37. The top curve 
is the Fourier spectrum of the pendulum twist from an 
Experiment I data run and shows the driven harmonic 
peaks as well as the free torsional oscillation (these data 
were not passed through the torsion filter). The effect of 
the pendulum as a low-pass filter is seen by the drop in 
the noise floor above the resonance at 17w. The smooth 
curve is the predicted thermal noise from damping inter- 
nal to the torsion fiber. The spectrum follows the curve 
well, indicating that we were operating near the ther- 
mal limit. It is interesting to note that velocity damping 
(from residual gas, for example) is negligible as it would 
produce a constant, as opposed to 1/f, background at 
low frequencies. The bottom curve shows that the angle 
readout noise was well below the "torque noise" associ- 
ated with actual movement of the pendulum. 
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2 5 10 20 50 

Attractor Rotation Harmonic 

FIG. 37: Fourier spectrum of the pendulum twist from an 
Experiment I data run. The horizontal axis is displayed in 
units of attractor rotation harmonic. The IOoj, 20a>, and 30u> 
signals are easily recognizable, along with the free torsional 
motion at 11 lo. The smooth curve shows the predicted ther- 
mal noise plus a white noise floor that arose from ADC "bit 
noise". The bottom curve, which shows the readout noise 
in the twist signal, is the spectral density recorded with the 
pendulum held stationary. The "bit noise" in this spectrum 
is negligible because a higher gain was used. 



APPENDIX B: CALCULATING NEWTONIAN, 
YUKAWA, POWER-LAW AND MASSIVE 
PSEUDOSCALAR FORCES BETWEEN 
CYLINDERS 

1. Approach 

We assume the test bodies are perfect cylinders with 
symmetry axes along z, radii and thicknesses a\, h\ 
and a 2 , h 2 , and uniform negative densities p\ and p 2 . 
Without loss of generality, we can take one cylinder to 
be centered at the origin and the other at coordinates 
(i, 0, s + [hi + h 2 }/2) as shown in Fig. 38. 

The Newtonian, Yukawa, power-law and massive 
pseudo-scalar force calculations (see Eq. 56 each required 
evaluating 6 integrals. We evaluated analytically as many 
of the integrals as possible, and did the remaining ones 
numerically. We expedited the calculations by making 
databases of the various forces evaluated on a grid of 
horizontal and vertical separations t and s (for Yukawa 
and massive pseudo-scalar forces the ranges A were an ad- 
ditional parameter) . The force between a pendulum and 
upper-attractor hole pair were calculated for 500 (100 for 
the non-Newtonian forces) values of t spaced quadrati- 
cally in the region 0.0 mm < t < 85 mm (this placed more 
points where Ft is highly peaked) for each of 100 s values 
in the relevant regime. The force from a lowcr-attractor 
hole was calculated for 0.05 mm + h < s < 11 mm + h 
where h is the thickness of the upper disk). The hor- 
izontal forces at arbitrary t and s were found by cubic 
spline interpolation of the database values as described 




FIG. 38: Coordinate system used in the force calculations, 
in Ref. [37]. 

In all cases, the x and z integrals were evaluated first 
to give the force per unit thickness between two infinitely 
thin rectangular sheets lying parallel to the x — z plane, 
as illustrated in Fig. 38. Because the force depends 
only upon the coordinate differences (x± — x<i, etc.), we 
transformed coordinates to x = x\ — X2, X = xi + X2, 
z = Z\ — Z2, and Z = z\ + z 2 . Then the x integrations 
gave 



——dx 1 dx 2 = 
ox Jf 



Vdx 



Vdx 



(Bl) 



where V is one of the potentials in the square brackets 
of Eqns. 56 and 

(B2) 



to 


= t- 


(4 


-vl) 1/2 - 


(0? 


-vi) 1/2 


tl 


= t- 


(4 


~vl) 1/2 + 


(a? 


-yl) 1/2 




= i + 


(4 


-vl) l/2 + 


(a? 


-yf) 1/2 


*3 


= i + 


(4 


-vl) 112 - 


(4 


-y\f 12 



Similarly, the z integrations gave 
Ixz = I dz 1 dz 2 I x 



(B3) 



(z-s Q )I x dz+ / (s 1 -s )I x dz+ / (s 2 -z)I x dz 



where sq — s, s\ — s+hi, s 2 = s+hi+h 2l and S3 = s+h 2 . 
Finally, the y integrals were evaluated numerically using 
a standard Romberg integration scheme [37]. Reflection 
symmetry about the x — z plane allowed us to integrate 
one half of one cylinder and double the result, giving 



F t (t, s) 



2a\a 2 



ir/2 j.n/2 
Tt/2 JO 



Gp\p2 

where we made the substitution yi — > otj sin 9 



I xz cos 0i cos9 2 d6i d9 2 . 

(B4) 
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2. Newtonian Force Calculation 



and p = y + z , we have 



All four Newtonian x and z integrals were evaluated 
analytically, yielding 



3 3 



4(t.j.') = -iEE(-i) <+j 



(B5) 



i=0 j=0 

i i ti / 2 2 \ l r<K r ij 

** r y + T7r(y s i) ln 



2t 4 s, In I 1 + ) + 2ys, tan" 1 ( 



where y = y\ — y 2 and r 2 - = i 2 + y 2 + s 2 . 



3. Yukawa Force Calculation 



F (2 

d 2 



13 



(3 

y 

G (3 

G 



k, I r 
2 ln 



+ pj tan 1 



- tan" 1 - 
P \P 



hi 



U 



p 2 VW- 
1 

- — tan 

U 



Pj 



G?> = 



1 



tan 



3p 2 r y 



(B9) 



We integrated Eq. B3 for the Yukawa interaction by 
parts and made the substitution p = 1/z. This left the 
following x and z integrals, along with the y integrals, to 
be evaluated numerically 



JTz(*,2/,s,A) = A 



*2 ftl 

dx- dx (B6) 

^3 ^ *0 

l/s /-l/ s 3 rl/s 

s / + (si-s ) / Hdp + s 2 I Hdp 

Ji/ ai 



l/si 



1/S3 



5. Massive Pseudoscalar Exchange Calculation 

In this case, only the integrations leading to Eqs. Bl 
and B3 could be done analytically. The 4 integrations 
needed to evaluate Eqs. Bl, B3 and B4 were done nu- 
merically. 



where 



H 



cxp • 



-V^ 2 + y 2 + i/p 2 



(B7) 



4. Power-law Force Calculation 



The power law potentials, V oc 1/r for k = 2 — 5, 
required three numerical integrations. In this case, I xz 
took the form 



/£ (t, y, *, A) ^g" 1 



i=0 



(fe) 



(B8) 



/SO /*S3 /■ 

G\ k) dz + ( Sl -s ) / G| fc) dB + s 2 / 
1 J Si J 8, 



Gf'dz + s 2 I G\ k) dz 

S3 



Letting y = y 1 - y 2 , r 2 = i 2 + y 2 + s 2 p 2 = y 2 + s 2 , 
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